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Confinement—deconfinement phase transition 
and gauge-invariant gluonic mass in Yang-Mills theory 

Kei-Ichi KondcQ 

Department of Physics, Graduate School of Science, Chiba University, Chiba 263-8522, Japan 

We give an analytical derivation of the confinement/deconfinement phase transition at finite 
temperature in the SU{N) Yang-Mills theory in the D-dimensional space time for D > 2. We 
elucidate what is the mechanism for quark confinement and deconfinement at finite temperature 
and why the phase transition occurs at a certain temperature. For this purpose, we use a novel 
reformulation of the Yang-Mills theory which allows the gauge-invariant gluonic mass term and 
calculate analytically the effective potential of the Polyakov loop average concretely for the SU(2) 
and SU{3) Yang-Mills theories by including the gauge-invariant dynamical gluonic mass. For D = 4, 
we give an estimate on the transition temperature Td as the ratio to the gauge-invariant gluonic 
mass M which has been measured on the lattice at zero temperature and is calculable also at finite 
temperature. We show that the order of the phase transition at Td is the second order for SU (2) and 
weakly first order for SU{3) Yang-Mills theory. These initial results are obtained easily based on 
the analytical calculations of the “one-loop type” in the first approximation. Then these results are 
identified with the initial condition in solving the flow equation of the Wetterich type to improve the 
initial results in a systematic way in the framework of the functional renormalization group. But the 
improvements do not change the initial results in an essential way except for some thermodynamic 
observables. We argue how the artifacts in the first approximation are eliminated to obtain the 
correct behaviors for such thermodynamic observables. 

PACS numbers: 12.38.Aw, 21.65.Qr 


I. INTRODUCTION AND MAIN RESULTS 

Quark confinement and chiral-symmetry breaking are 
the main subjects to be investigated for understanding 
the various phases in the gauge theory for strong in¬ 
teractions, namely QCD at finite temperature and den¬ 
sity, see e.g., [l| for recent texts and 0 for recent re¬ 
views. In a previous paper Q, we have proposed a 
theoretical framework to obtain a low-energy effective 
theory of QCD towards a first-principle derivation of 
confinement/deconfinement and chiral-symmetry break¬ 
ing/restoration crossover transitions at finite tempera¬ 
ture. The basic ingredients are a novel reformulation 
[5 of Yang-Mills theory and QCD based on new vari¬ 
ables originated fro m [5| -[l4|. and the flow equation of 
the Wetterich type in the framework of the func¬ 
tional renormalization group (FRG) [l^ as a realization 
of the Wilsonian renormalization group [T^ . In fact, we 
have demonstrated that an effective theory obtained in 
this framework enables us to treat both transitions si¬ 
multaneously on equal footing from QCD. The resulting 
effective theory in simple but non-trivial approximations 
is regarded as a modified and improved version of nonlo¬ 
cal Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) 
models proposed in [Tsl - I^ extending the original (local) 
PNJL model [2l|. See e.g., for later developments 

of the nonlocal PNJL models. 

A novel feature of the resulting effective theory is that 
the nonlocal Nambu-Jona-Lasinio (NJL) coupling de- 
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pends explicitly on the temperature and Polyakov loop, 
which affects the entanglement between confinement and 
chiral symmetry breaking, in addition to the well-known 
coupling in the conventional PNJL model between the 
Polyakov loop and quarks coming from the cross term 
generated through the covariant derivative in the quark 
sector. The chiral symmetry breaking/restoration transi¬ 
tion is mainly controlled by the nonlocal NJL interaction. 

On the other hand, the confinement/deconfinement 
transition in the pure gluon sector [2^ is described by 
the nonperturbative effective potential for the Polyakov 
loop average which is obtained in a nonperturbative 
way put forward by [1^, [2^ in the framework of FRG 
(See also [jol - l^ i. At present, however, the FRG studies 
of the Yang-Mills theory and QCD rely heavily on hard 
numerical works and the outcome is obtained only in the 
numerical way. This fact unables everyone to reproduce 
the FRG results and to understand the results in a physi¬ 
cally transparent manner. Therefore, a simple analytical 
derivation of the results is desired to understand such 
nonperturbative phenomena from the first principle. 

In this paper, we focus on the confine¬ 
ment/deconfinement phase transition in the pure 
gluon sector for SU{N) Yang-Mills theory. We 
demonstrate that the essential features on the confine¬ 
ment/deconfinement phase transition summarized below 
can be obtained in a simple analytical way without 
hard numerical works, once we take into account a 
gauge-invariant and dynamical gluonic mass M which 
is allowed to introduce in the reformulation of the 
Yang-Mills theory. In fact, we have already emphasized 
the importance of such a gluonic mass in the previous 
papers [J, but have not exhausted the outcome yet. 
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The results obtained in this paper will be applied to the 
chiral-symmetry breaking/restoration and its crossover 
to confinement/deconfinement in QCD in a subsequent 
paper. 

For this purpose, we use the reformulation of the Yang- 
Mills theory which allows one to introduce a gauge- 
invariant “mass term” for a specific gluonic degree of free¬ 
dom called the remaining field see e.g., 0 for a 

review. Such a gluonic mass has already played the very 
important role in quark confinement at zero temperature 
to explain and understand the “Abelian dominance” in 
the Maximally Abelian gauge gauge^ which is replace d by 
the gauge-independent restricted field dominance [42 h 4^ 
in our terminology [l^ . 

First of all, notice that the mechanism of the dynami¬ 
cal mass generation for the gluon field has been already 
proposed and that the dynamical gluonic mass generation 
has been shown to occur at zero temperature in [T^ , 

see also for the related works. 

0. The gauge-invariant mass M for the remaining field 
can be generated dynamically through the 
gauge-invariant vacuum condensation of mass di¬ 
mension two:^ 

$ := (= (2tr[trp.r^]), (i.i) 

which occurs due to the quartic self-interactions 
among the gluons represented by the remaining 
fields in the Yang-Mills theory. The dynamical glu¬ 
onic mass M is obtained from the minimum of the 
effective potential of the vacuum conden¬ 

sate $, which is also written as 14ff(M). Another 
way of understanding the mass term is also given 
from the viewpoint of the gluonic Higgs field, which 
can be elucidated only in our formulation. See sec¬ 
tion ESI 


finite temperature. At the same time, we want to calcu¬ 
late the Polyakov loop average L to discuss the confine¬ 
ment/deconfinement transition. Therefore, we need to 
calculate the simultaneous effective potential 
as a function of the two variables M and L. It is 
confirmed by comparing and [13 that the gauge- 
invariant vacuum condensation of mass dimension two 
$ = can be related to the well-known gauge- 

invariant gluon condensation of mass dimension four, i.e., 
responsible for the trace anomaly, which de¬ 
termines the non-perturbative vacuum, see also In 
this paper, however, we treat the mass M just as a con¬ 
stant without the temperature dependence by restricting 
to the effective potential t4ff(L) of the Polyakov loop av¬ 
erage L alone for simplicity. Hence, M is equal to the 
value at zero temperature. The result of the effective po¬ 
tential I4ff($, L) will be reported in a subsequent work. 

The following results are obtained based on an ana¬ 
lytical calculation of the effective potential Ves{L) of the 
Polyakov loop average L alone in the SU (2) and SU (3) 
Yang-Mills theories at finite temperature T in I? = 4 di¬ 
mensions by including the gauge-invariant and dynamical 
“gluonic mass” M. 


1. There exists a confinement/deconfinement phase 
transition at a critical temperature in the re¬ 
spective Yang-Mills theory at finite temperature 
T signaled by the Polyakov loop average {L{x)), 
i.e., non-vanishing {L{x)) ^ 0 for high tempera¬ 
ture T > Td^ and vanishing {L{x)) = 0 for low 
temperature T < Td- ^ The Z{N) center symmetry 
which is spontaneously broken at high temperature 
restores at low temperature.^ 

2. The critical temperature Td is estimated in the form 
of the ratio to the dynamical gluonic mass M in the 
respective Yang-Mills theory: 


The above ideas enable us to calculate the gauge- 
invariant dynamical gluonic mass M also at finite tem¬ 
perature in our reformulation. In fact, the temperature 
dependence of the dynamical mass M{T) is obtained 
from the minimum of the effective potential Ves{M) at 


Td/M =0.34 for S't/(2), 

Td/M =0.36 for S't/(3). (1.2) 

It should be emphasized that this ratio is gauge- 
independent. To obtain the critical temperature 
Td, we need to know the value M of the gluonic 


^ The Abelian dominance in the low energy regime in the Yang- 
Mills theory was first proposed in l35l and was confirmed for the 
string tension by the numerical simulations on the lattice [3^ in 
the Maximally Abelian gauge [371 . which is a realization of the 
Abelian projection l34l . The magnetic monopole dominance was 
also confirmed for the string tension by the numerical simulations 
on the lattice in the Maximally Abelian gauge l38l . The Abelian 
dominance was also shown for the gluon field propagators l4(j . 
|41| . These results suggest the dual superconductor hypothesis 
as a promising mechanism for quark confinement l39l . 

^ The condensate 4’ is a gauge-invariant version (which is made 
possible in our formulation) of the BRST-invariant vacuum con¬ 
densation of mass dimension-two obtained from the on-shell 
BRST invariant operator of mass dimension two proposed in |49|| . 


® On a lattice, a rigorous proof for the existence of SU(N) and 
U (N) gauge theory in d > 3 dimensions was given by l26l . 

The center symmetry corresponds to the aperiodic gauge trans¬ 
formation: U{t + T~^,x) = U{t,x)g (U,g G G) such that the 
gauge transformed field is periodic: {t + T' ^ (i , x) 

(the periodicity of the gauge field is preserved under the gauge 
transformation + T“^,£c) = ®)) f53| . Under such an 

aperiodic gauge transformation, the action is invariant, while the 
Polyakov loop operator is not invariant: L{x) gL[x). Such a 
set of elements g constitutes a discrete subgroup of G called the 
center. For G = SU{N), g = zl must satisfy gg'^ = 1 and 
det g = 1, which yields 22 * = 1 and = 1. From this observa¬ 
tion, we find Center(SU(N)) = Z(N). 
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mass.® The values of the gluonic mass M have 
been measured on the lattice at zero temperature 
T = 0 by Shibata et al. 

M{T = 0) =1.1 GeV for SU(2), 

M{T = 0) =0.8 - 1.0 GeV for SU{3). (1.3) 

A naive use of these values of M leads to the esti¬ 
mate on T^-. 

Td =374 MeV for SU(2), 

Td =288 - 360 MeV for 5't/(3). (1.4) 

Incidentally, the numerical simulations on a lattice 
give the values [H^: 

Td =295 MeV for 5t7(2), 

Td =270 MeV for 5t7(3), (1.5) 

while the continuum a ppr oach, e.g., the most recent 
FRG studies give [1^ [sH 

Td =230 MeV for S'17(2), 

Td =275 MeV for S'17(3). (1.6) 

3. The order of the phase transition at Td is the sec¬ 
ond order for SU{2) and (weakly) first order for 
SU{3) Yang-Mills theory. This result is shown to 
be consistent with the standard argument based 
on the Landau theory of phase transition using the 
expansion of the effective potential 14ff(L) into the 
power series of the Polyakov loop average L as the 
order parameter. In particular, the first order tran¬ 
sition in the SU (3) Yang-Mills theory is induced by 
the cubic term L® of the Polyakov loop average L 
in the effective potential Vee{L). 

4. The mechanism for quark confinement or decon¬ 
finement at finite temperature is elucidated with¬ 
out detailed numerical analysis in this framework 
by taking into account the gluonic mass M. In high 
temperature T ^ M the gluonic mass M becomes 
negligible and all the relevant degrees of freedom 
behave as massless modes, and the effective poten¬ 
tial can be calculated in the perturbation theory so 
that the minimum of the effective potential Ves{L) 


® Our estimate on is indeed a little bit higher than expected at 
present. But this is based on the value of the mass M obtained at 
zero temperature T = 0. The gluonic mass M should depend on 
the temperature T. The mass M should be determined in a self- 
consistent way, not just a given parameter. Indeed, if the mass 
M decreases as the temperature increases: M{T > 0) < M{T = 
0), then the initial value reproduces a better result than the 
naive estimate. Therefore, our approach has the potential to give 
better numerical estimate on without further improvements. 
The direct measurement of the gluonic mass M on the lattice at 
finite temperature is under way 1531 . 


is given at the non-vanishing Polyakov loop average 
L ^ 0 implying deconfinement [H, . Whereas 

in low temperature T M the “massive” spin- 
one gluonic degrees of freedom (i.e., two transverse 
modes and one longitudinal mode) are surpressed 
and the remaining unphysical massless degrees of 
freedom (i.e., a scalar mode, and ghost-antighost 
modes) become dominant. Consequently, the sig¬ 
nature of the effective potential Ves{L) is reversed 
so that the minimum of the effective potential is 
given at the vanishing Polyakov loop average L = 0 
implying confinement.® 

5. The above results are shown using the first approx¬ 
imation based on the analytical calculations of the 
“one-loop type” (which is different from the one- 
loop calculation in perturbation theory). This re¬ 
sults of the first approximation offer an effective 
starting point for the more systematic analysis of 
the non-perturbative studies.^ These initial re¬ 
sults are regarded as the initial condition in solving 
the flow equation of the Wetterich type and they 
can be improved in a systematic way in the FRG 
framework accordir^ to the prescription given in 
the previous paper |3| where the crossover between 
confinement/deconfinement and chiral symmetry 
breaking/restoration has been discussed from the 
first principle, i.e., QCD, without explicitly intro¬ 
ducing the gluonic mass. [ Remember that the first 
approximate solution of the Wetterich equation is 
given by the “one-loop type” expression with the 
additional infrared regulator term which plays the 
role of the mass term in a certain sense. ] But, the 
FRG improvement does not change the above con¬ 
clusions in an essential manner. The above Td gives 
a lower bound on the true critical temperature Tc, 
since the flow evolves towards enhancing the con¬ 
finement, under the assumption that M does not 
change so much along the flow. 

6. Remark: We must be cautious in treating the ther¬ 

modynamic observables, which needs the value of 
the absolute minimum 14 = t4ff(Tmin) of the ef¬ 

fective potential 14ff, he., the vacuum energy. We 
do not need such information to derive the above 
results which are obtained only from the location 


° This observation is in line with the general arguments given in 
in the FRG and agrees with the statement given in [s^ . Is^l 
in the approach [s^ . 

^ The first approximation does not depend on the gauge coupling 
constant explicitly, although the coupling constant dependence 
could be included into the result through the gluonic mass implic¬ 
itly. The explicit dependence appears if we calculate the gluonic 
mass from the simultaneous effective potential The 

improvement by the FRG approach depends explicitly on the 
coupling constant. 
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^min of L giving the minimum 

dV,s{L) 


Kff (^11 


dL 


= 0 . 


L=L„ 


(1.7) 


The thermodynamic pressure P{T) = —= 
— t4ff (Lmin(7’)) remains positive in the low- 
temperature confined phase L = 0 in the first ap¬ 
proximation of our formulation, in sharp contrast 
to the positivity violation reported in the preceding 
work at one loop [H, . For the entropy density 

S{T) := we find the positivity violation near 

the critical temperature and need the improvement 
of the naive first approximation. We discuss the 
theoretical and physical reasons for these artifacts 
in Section flV PI 


approximation. Moreover, we discuss how the initial re¬ 
sults obtained in the first approximation are improved 
by using the flow equation of the Wetterich type in the 
framework of FRG in a systematic way. Finally, we point 
out the artifacts in treating the thermodynamic quanti¬ 
ties, e.g., the pressure and the entropy, and propose the 
possible resolution in our framework. 

In Section V, we show the existence of the conflne- 
ment/deconfinement phase transition at finite tempera¬ 
ture in SU (3) Yang-Mills theory. We show that the criti¬ 
cal temperature Td is estimated as the ratio Td/M to the 
dynamical gluonic mass M, and that the phase transition 
is discontinuous, namely, the weakly first order. 

The final section is devoted to conclusions and dis¬ 
cussions. Some technical materials are collected in the 
Appendices. 


Besides the numerical simulations on the lattice, there 
are other approaches, see e .g., Ibp-lSSj. Among them, 
especially, the authors of [^, have introduced a 
different kind of gluonic mass term in the gauge-fixed 
Yang-Mills theory at finite temperature and have in¬ 
vestigated the effect of the mass term on conflne- 
ment/deconfinement phase transition. They have found 
that the phase transition is quite well described by the 
one-loop calculations in the perturbation theory, once 
the gluonic mass is introduced to the Yang-Mills theory. 
Their works are very interesting in its own right, but 
quite surprising. One must answer what is the meaning 
of the gluonic mass and why the one-loop calculation is 
enough (although they tried to improve the one-loop re¬ 
sult by including the two-loop result M)- We will give a 
partial answer to these questions from our point of view. 
It should be remarked that their mass term is somewhat 
similar to ours at first glance, but its theoretical origin 
and the content are totally different from ours. 

This paper is organized as follows. 

In Section II, we give the reformulation of the SU(N) 
Yang-Mills theory written in terms of the new variables. 
We show that the gluonic mass term is introduced in 
the gauge invariant way in the reformulated Yang-Mills 
theory. We discuss also the meaning of the gluonic mass 
term. 

In Section III, we discuss how the Polyakov loop oper¬ 
ator is expressed in terms of the new variable. Then we 
explain our standpoint to give a prescription of calculat¬ 
ing the effective potential of the Polyakov loop average 
in the reformulated SU(N) Yang-Mills theory. 

In Section IV, we show the existence of conflne- 
ment/deconfinement phase transition at finite tempera¬ 
ture in SU{2) Yang-Mills theory by examining the effec¬ 
tive potential of the Polyakov loop average. This result 
clarifies the mechanism of the phase transition. We show 
that the critical temperature Td is estimated as the ra¬ 
tio Td/M to the dynamical gluonic mass M, and that 
the phase transition is continuous, namely, the second 
order. These are the initial results due to the analytical 
calculation of the “one-loop type” obtained in the first 


II. REFORMULATING YANG-MILLS THEORY 
USING NEW VARIABLES 

A. Gauge-covariant decomposition of Yang-Mills 
field 

In this section, we give a brief introduction to the refor¬ 
mulation of the Yang-Mills theory needed in what follows, 
see e.g., a review [J for the details. In this paper we con¬ 
sider the SU{N) Yang-Mills theory in the H-dimensional 
space-time. We use the gauge-covariant decomposition of 
the SU{N) Yang-Mills field into two pieces f^(cc) 

and {p, = 0,1,..., D — 1): 

s^^{x) = y^{x) + 3S^{x) e su{N) := Lie{SU{N)), 

(11.1) 

where := He{G) denotes the Lie algebra of a Lie group 
G. When the original Yang-Mills field s^fi{x) obeys the 
ordinary gauge transformations given by 

s^^{x) £//,{x) := Uix)[j!/f,ix)+ig-^d^]Uix)-\ 

(11.2) 

the first piece Tfiix) called the restricted field and the 
second piece 3S^(x) called the remaining field are re¬ 
quired to obey the gauge transformation: 

y^{x) T/ix) := Uix)['r^{x)+ig-^df,]U{x)-\ 
^^{x) 3T/{x) := U{x)3T^,{x)U{x)-^. (11.3) 

Therefore, we have the same form of the decomposition 
after the gauge transformation: 

= r/ix) + s:/{x) G su(N). (11.4) 

Such a decomposition can be constructed by introduc¬ 
ing a Lie algebra valued field nix) called the color di¬ 
rection field or color field for short which is supposed 
to obey the gauge transformation in the adjoint repre¬ 
sentation: 

n{x) n'{x) := U{x)n{x)U{x)~^ G Lie{G/H). (11.5) 
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Here H is a subgroup of G called the maximal stabil¬ 
ity subgroup There are a number of options 

depending on the choice of H. For G = SU{N), the max¬ 
imal stability subgroup H is equal to the maximal torus 
subgroup H = H ■= in the m^Lximal option, 

and H = U{N — 1) in the minimal option. 

In this paper we discuss only the maximal option of 
the SU(N) Yang-Mills theory and omit other options, see 
[1) Hi- The group G = SU(N) has the rank r = N — 1. 
In the maximal option, it is possible to construct a set of 
r Lie algebra 1^-valued fields nj{x) (j = 1, • • • , r) by the 
repeated multiplication of the original color field n{x)'. 

nj{x) =nf{x)TA& Lie{G/H) (j = l,..-,r), (II.6) 

where Ta {A = I,..., dimG = — 1) are the generators 

of su{N): Ta = ^(Ja with a a being the Pauli matri¬ 
ces for G = SU{2) and Ta = with being the 
Gell-Mann matrices for G = SU{3). The color fields are 
orthonormal; 

nj{x) ■ nk{x) := 2tr{nj{x)nk{x)) = nf{x)nk{x) = 6jk, 
j,fce{I,2,--- ,r}, (II.7) 


The commutativity (III.8|) is also satisfied, since Hj are 
the Cartan subalgebra obeying 

[H„Hk]=0, j,k€{l,2,--- ,r}. (11.14) 

Once a set of color fields nj{x) satisfying the above 
properties is given, the respective pieces '7^(x) and 
^fj,{x) of the decomposition are uniquely determined by 
imposing the conditions called the defining equation: 

(I) all color fields rij (x) are covariantly constant in the 
restricted background field 'f'fj,{x): 

0 = ^f,['V]nj{x) =df,nj(x) - ig[y^{x),nj{x)] 

(j = l,2,--- ,r), (11.15) 

(II) the remaining field ^fi{x) is orthogonal to all 
nj{x)-. 

0 = ST^{x) ■ nj{x) :=2tT{^fj,{x)nj{x)) = S^^{x)nf{x) 
(j = l,2,.-- ,r). (11.16) 

The defining equation (I) follows from 
(F) The single color field n{x) is covariantly constant in 
the background field ^^(a:): 


and they mutually commute: 


0 = ^fj,[y]n{x) = d^n{x) - ig[yfj,{x),n{x)\. (11.17) 


[nj{x),nk{x)\=Q, j, fc G {1, 2, • • • , r}. (II.8) 

Therefore, all rij (x) have the same gauge transformation 
as n(x): 

nj(x) —>■ n'^{x) := U{x)nj{x)U{x)~^ 

G Lte{SU{N)/U{l)^-^) = 

(II.9) 

For G = SU{3), we introduce the two color fields denoted 
by ria and ng: 

n 3 (a:) = n^{x)TA, ns{x) = n^{x)TA G Lie{SU{3)/U{Vf). 

(II.IO) 


For G = SU{3), n is constructed as a linear combination 
of ng and ng. A simple choice is n{x) = n^{x). Then 
ng is constructed from n^. Indeed, the two color fields 
are related as 


n3(a:)n3(a:) = 


(II.ll) 


Such color fields nj{x) are constructed using the ad¬ 
joint orbit representation from the generators Hj of the 
Cartan subalgebra of 1^: 

nj{x) = U\x)HjU{x) G Lie{G/H), j G {1, 2, • • • , r}. 

( 11 . 12 ) 


In fact, the fields nj(x) defined in this way satisfy the 
ortho-normality condition (III.71) . since 


The defining equation (II) is also given by 

(IF) ^^{x) does not have the iF-commutative part 

^^{x)h '■= ^^{x) — [nj{x), [nj{x), ^^{x)]]: 

0 = 3T^{x)h ^^{x) = [nj{x), [nj{x), ^^(x)]]. 

(11.18) 

In what follows, the summation over the index j = 1,..., r 
should be understood when it is repeated, unless other¬ 
wise stated. It is possible to show [J the equivalence 
between (111.161) and (111.181) . 

By solving the defining equations, the decomposed 
fields Y(i(x) and STfj_{x) are determined uniquely: 

= - ig~^[nj{x), %[£/]nj{x)] G Lie{G/H), 

Tk-ix) =‘^/x(a;) + -^kix) G Lie{G), 

^fj.{x) =nj(x)(nj(x) ■ s^k^x)) = nj{x)‘^i{{x) G Lie{H), 
^kix) G Lie{G/H). (11.19) 

Thus, once a full set of color fields nj(x) is given, the 
original gauge field has the unique decomposition called 
the Cho-Faddeev-Niemi decomposition H. [l3|. 
In this stage, the decomposed fields are written in terms 
of nj{x) and £^^{x). For SU{2), in particular, we have a 
single color field n(x) and the decomposition reads 

■^ki^) = - WymIMx), ^ Lie{SU{2)/U{l)), 

Tkix) =‘^ki^) +^ki^) ^ Lie{SU{2)), 

‘^fj,{x) =n{x){n{x) ■ ^;i(x)) = n(x)'^^(x) G Lie{U{l)), 
^kix) =ig~^[n{x),d^,n{x)] G Lie{SU{2)/U{l)). 

( 11 . 20 ) 


nj{x) ■ nk{x) = 2tr{HjHk) = Hj ■ Hk = Sjk. (11.13) 
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This is the Cho-Duan-Ge decomposition mm- 

M- 

The advantages of the decomposition are as follows. 

(a) [restricted field dominance] The original Wilson loop 

operator and the Polyakov l oop operator are exactly re¬ 
produced from alone [gl, 113, : 

WcH = VPcm, L4^] = L4i^], (11.21) 

(b) [gauge-invariant field strength] The gauge-invariant 

field strength is obtained from the field strength of 
the restricted field \y\ := [y^, 

in the rij direction [J: 

=ti{nj{x)^f,^['r]{x)}. ( 11 . 22 ) 


B. Reformulation of Yang-Mills theory 

The goal of the reformulation is to change the original 
field variables into the new field variables.® For 

this purpose, the color field n{x) G Lie{G/H) must be 
written in terms of the original £/^(x) G Lie{G). This is 
achieved by solving the reduction condition x(a:) = 0 
for a given sy^{x). A choice of the reduction condition is 

•= [nj{x),^^[sy]2!fj,[sy]nj{x)\ G Lie{G/H). 

(11.23) 

Thus, all the new field variables have been written in 
terms of the original variables 

In the original Yang-Mill theory, the average of F is 
written (omitting the gauge fixing of G) as 

{F[£/])ym = Zyl J (11.24) 


® The reduction condition is necessary for the reformulated gauge 
theory written in terms of the new field variables based on change 
of field variables to be equivalent to the original gauge theory. It 
has been shown El, El that the reduction condition is regarded 
as the gauge-fixing condition of breaking the enlarged gauge sym¬ 
metry G X G/H into the original gauge symmetry G. Even after 
imposing the reduction condition, therefore, the original gauge 
symmetry is retained in the reformulated gauge theory. The 
reduction condition is given by the global minimum condition 
of a functional defined as the integration over the whole space- 
time, which is minimized under the gauge transformation of the 
enlarged gauge symmetry El El. Therefore, the reduction con¬ 
dition is free from the Gribov copy problem in principle. In fact, 
we can transfer this framework to the lattice to perform the full 
non-perturbative studies The reduction condition and 

the associated “Faddeev-Po^v” determinant (which was called 
the Shabanov determinant Q) given in Section II are just the 
local representations corresponding to the local minimum given 
in the framework of the local field theory. For the details, see the 
original references EMI or the review [J|. The gauge fixing of 
the original gauge symmetry G is discussed in the next section 
III. 


In the reformulated Yang-Mills theory, the average of F 
is rewritten as 

(F’K])ym' =z-l,, J 

X jrJ’]. (11.25) 

(See [liHil for SU{2), M for SU{N), which correspond 
to section 4.4 for SU{2), and section 5.6 for SU{N) in 
the review Q.) Here (i) X = 0 is the reduction condition 
rewritten in terms of the new variables: 

X:=x[n,^,^]:=D>^[y]^^. (11.26) 

This constraint can be incorporated into the Lagrangian 
by introducing the Lagrange multiplier field, namely, the 
Nakanishi-Lautrup field ^{x): 

5(x) = J (11.27) 

= 2xv[jrs>^[r]SF^]. (11.28) 

(ii) App is the “Faddeev-Popov” determinant associated 
with the reduction condition. The precise form is ob¬ 
tained by the BRST method as 

Ag*^[n,'r, JT] = det{-^^[r-jr]^'^[r-Hjr]}. (11.29) 

(See 0, section 4.6 and Appendix E for SU{2), and 
section 5.8 and Appendix H for SU{N) in Q.) The de¬ 
terminant can be incorporated into the Lagrangian by 
introducing the “ghost” ri{x) and “antighost” held fj{x) 
as 

Ap*^[n,'r, ^] = j 

jSfpp = ifi^{-9Fy - S:]9^[r + (11.30) 

(iii) J is the Jacobian associated with the change of vari¬ 
ables. By a suitable choice of the basis in the color space, 

J = l, (11.31) 

irrespective of the choice of the reduction condition. (See 
0 for SU{2), 0 for SU{N), which correspond to sec¬ 
tion 4.5 for SU{2), and section 5.7 for SU{N) in the 
review Q.) 

Notice that the reformulated Yang-Mills theory retains 
the gauge symmetry of the gauge group G = SU{N) 
of the original Yang-Mills theory at any stage even af¬ 
ter imposing the reduction condition and introducing the 
“Faddeev-Popov” determinant. 



original —>■ reformulated 

held variables 


action 

-a SYyiin,^^, 

measure 


operator 

FH ^ jrj>] 


Remarkably, this reformulation allows us to introduce 


the gauge-invariant “mass term” for the remaining 
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field which is a piece obtained by the gauge-covariant 
decomposition -p of the original gauge field 

(11.32) 

In fact, the numerical simulations on the lattice exhibit 
the dynamical mass generation for the remaining 
field at zero temperature, see [i^ for SU{2), [i^ for 
51/(3) maximal option, and [i^ for SU{3) minimal op¬ 
tion: 

M{T = 0) ~ 1.1 GeV for SU(2), 

M{T = 0) ~ 0.80 ^1.0 GeV for SU{3) maximal option, 

M{T = 0) ~ 1.1 ^ 1.15 GeV for SU{3) minimal option. 

(11.33) 

In the reformulated Yang-Mills theory, thus, the aver¬ 
age is obtained as 

(FK])ym' =^ym' / 

(11.34) 

where 

J/, ?7, ??], 

V, V] =-^YM[n, 

-|-.ifred[w,'^, 

+ (11.35) 

In the reformulated Yang-Mills theory, the Polyakov loop 

average is obtained as 

{Lx[£/])ym' 

=^YM' J 

X J (11.36) 

We can show that the Lagrangian density of the 
SU{N) Yang-Mills theory 

^YM = -itr(^^,.K]^'^"H), (11.37) 

is decomposed into the form: 

a:''^ 

+ o{a:^), 

+ (11.38) 


where is the non-Abelian field strength of the 

restricted field In the Yang-Mills Lagrangian, the 
terms linear in the remaining field do not appear. In 
what follows, we neglect all terms cubic and quartic in 
the remaining field a^p,. In particular, the cubic terms do 
not appear for the 51/(2) group and for the SU{3) group 
in the minimal option. 

C. A novel viewpoint for the gluonic mass term 

A novel viewpoint for understanding the gluonic mass 
term is as follows. For simplicity, we discuss only the 
51/(2) case. Then the gauge-invariant mass term (|II.32I) 
is rewritten in terms of the original variables s^p. 

S£^=M\x{{sap-'fpY} 

=MhT{{£/p - CpU - ig-^[dpn,n]f} 

(11-39) 

with the understanding that the color field n is expressed 
in terms of the original gauge field ji/p by solving the re¬ 
duction condition. Therefore, Yp (or Cp and n) plays the 
similar role to the Stuckelberg field to recover the local 
gauge symmetry. Note that Cp, n and a’p are treated as 
independent variables after the change of variables and 
the mass term is a polynomial in the new variable a'p, 
although they might be non-linear composite operators 
of the original variables £/p. 

We can identify the color field n{x) with the gluonic 
Higgs field (j>{x): 

(j){x) = Mn(x) G Lie{SU{2)/U{l)), (11.40) 

the mass term is regarded as the kinetic term for the 
non-linear sigma model: 

=g-MiD,Wr} + 42tr(<(.2) - M^} 

= ^5ym(-^a<H^) • {Dp[£/](j)) +u{(j3-(l)- M^), 

(11.41) 

where u{x) is the Lagrange multiplier field for incorpo¬ 
rating the constraint: 

(j)ix) ■ (j)ix) - = 0. (11.42) 

Alternatively, the mass term is regarded as the limit A —>■ 
oo of the model: 

■^m = ^'?ym(-^a*H<(') • + A((() • ()) - M2)2. 

(11.43) 

The dynamical generation of the gluonic mass is equiv¬ 
alent to the gluonic Higgs phase. Thus, the Yang-Mills 
theory with the “mass term” for the remaining field a^p 
is identified with the Yang-Mills-Higgs model with the 
gluonic Higgs field (j){x) in the Higgs phase 




The proposed mass term (IIL32I) for the gluon should 
be compared with the conventional gauge-invariant mass 
term of Kunimasa-Goto type [i^ : 

=Mhy{{^^ - ig-^Ud^U^f} 
=MhT{{UD^[£/]U'<f}, U{x) = 

(11.44) 

This mass term is non-polynomial in the Stiickelberg field 
X{x)- This fact makes the field theoretical treatment very 
difficult. 


D. Mechanism of dynamical gluonic mass 
generation 


Stokes theorem for the Wilson loop operator (III.211) and 
that the Wilson loop average is entirely estimated by the 
restricted field alone. 

The dynamical gluonic mass has various implications 
to the nonperturbative aspects of the Yang-Mills theory, 
e.g. . the glueball mass [^, the stability of the vacuum 
[4^ and to the low-energy effective theory for gluon con¬ 
finement Notice that there is an unsolved problem 
of a physical unitarity for the massive vector model . 

III. EFFECTIVE POTENTIAL OF THE 
POLYAKOV LOOP AVERAGE 

A. Polyakov loop operator 


It is possible to argue that there occurs a novel vac¬ 
uum condensation of mass dimension two for the 
field i.e.,® 


(11.45) 

and that the field acquires the mass dynamically 
through this condensation. This is a gauge-invariant ver¬ 
sion (which is made possible in our formulation) of the 
BRST-invariant vacuum condensation of mass dimension 
two obtained from the on-shell BRST invariant operator 
of mass dimension two proposed in (4^ . A naive way 
to see this is to apply the mean-field like argument or 
the Hartree-Fock approximation to the four-gluon inter¬ 
action, i.e., the quartic self-interaction among which 
leads to the gauge-invariant mass term for gluons and 
a gauge-invariant gluonic mass M : 


= - X ^.) • X a:'') 

=^M^a-^-a-f^, M^^‘^gl^{-a-fa-p^). (ii.46) 


The more detailed and systematic analyses have been 
given in [l^ . Consequently, the gauge-invariant mass 
for the “off-diagonal” gluon a’fj, is generated. Then the 
a"^ gluon modes decouple in the low-energy (or long¬ 
distance) region below the mass scale M. Consequently, 
the infrared “Abelian” dominance or restricted field 
dominance for the large Wilson loop average follows im¬ 
mediately from the fact that the Wilson loop operator is 
written in terms of alone according to the non-Abelian 


® We adopt the Minkowski metric = diag(l, —1, —1, —1). Af¬ 
ter the Wick rotation to the Euclidean region, the Minkowski 
metric tensor is replaced by = diag(—1, —1, —1, —1). 

Therefore, we have ——>■ > 0. 


The SU{N) Polyakov loop operator L^I-b/] is de¬ 
fined by taking the trace of the holonomy operator 
Ps^{x): 


iccK] := tr(P^(a;))/tr(l), 

a/T 

Psi^{x) = a e-xp ig / (irjYo(a;,T) 


e SU{N), 

(III.l) 


where a is the path ordering. The Polyakov loop op¬ 
erator Lx[.b/] is gauge invariant: Lx[-b/'] = Lx[-a]. In 
our reformulation, in can be replaced by the 

restricted field exactly: 


Lxi'P] := tr(Py(a;))/tr(l), 

nl/T 

P-r{x) = a exp 


ig / dT%{x,T) 
Jo 


e SU{N). 

(III.2) 


The Polyakov loop operator LxlP] is also gauge invari¬ 
ant: Lx[y] = Lx[^, and coincides exactly with the 
original one Lx [a'] [M . 1^ : 


Lx[£/] = Lx[r]. (III.3) 


B. Gauge fixing 

The reformulation given in the above is gauge invariant 
in the sense that the original gauge symmetry for the 
gauge group G = SU{N) in the Yang-Mills theory is 
retained at any stage. Now we proceed to discuss the 
gauge fixing in the reformulated Yang-Mills theory. 

The choice of the gauge does not change the aver¬ 
age value of the gauge invariant quantity such as the 
Polyakov loop average. In the actual calculations in this 
paper, therefore, we take a special gauge: the color field 
has the uniform direction of the Cartan subalgebra corre¬ 
sponding to the maximal torus subgroup PI = P(I)'^“^, 
which facilitates the calculation of the Polyakov loop 
average. This is equivalent to perform the local gauge 
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transformation to diagonalize the color field. Then the 
color fields nj{x) are no longer the local field variables 
and the integration over n in the path integral formula 
derived above becomes unnecessary. For SU{2), the re¬ 
sulting color field is chosen to be 

n'{x) = y ^ n^{x) = (III.4) 

Then the restricted field is given by 

=^'f,ix)n'{x)+ig-^[n\x),d^n\x)] 

=n(^)y- (ni.5) 

For SU (3), the color field is taken to be a linear combina¬ 
tion of the two diagonal generators Hi and H 2 belonging 
to the Cartan subalgebra: 

« 3 ( 2 ;) = y. = ^ ^ n'f{x) = 5f. (III. 6 ) 

Then the restricted field is given by 

In this way, the restricted field is chosen to be in the 
Cartan subalgebra. In what follows, we omit the prime. 
The restricted field is separated into the background part 
and the quantum fluctuation part: 

Y^{x)=f^{x) + %{x). (III.8) 

which is realized by separating (x) into the the back¬ 
ground part and the quantum fluctuation part. For later 
convenience, we choose the specific background = 

g~^Tip 6 ij,o of Y’ij,{x) for SU{2), and #^(x) = g~^T(fij 6 fj,o 
oiY’l'ix) for SU{3). 

For SU{2), 

Yfj.ix,T) = g~^T(pSfj,0Y 

^ r/(x,r) = g-^Tip5^oSA3 + (III.9) 

For 51/(3), 

Yf,ix,T) = g-^Tip3Sf,0Y + 9~^TipsSfj,QY + 

'^^ix,T) = g~^TSf,o{ip3SA3 + VsSas) + '^^(XjT), 

(III.IO) 


where ip, ips and ips are dimensionless quantities. 

We take the approximation in which the quantum fluc¬ 
tuation parts ly are neglected, ly <C 1 . Then the 
holonomy operator P{x) takes the simple form without 
the path ordering: 


P =exp 
P = exp 


. 0-3 


lip3- 


G [/(!) C 5t/(2) 
. A 


■ i(ps- 


G U{1) X C/(l) C 5t/(3), 

(III.ll) 


where <73 is the diagonal Pauli matrix and A 3 , Ag are the 
diagonal Gell-Mann matrices. The SU (2) Polyakov loop 
operator L becomes a real-valued function of an angle p: 


L{p) ■= = 2 


itr{exp[*^^]} 


= COS-G 


(III.12) 


and the SU{3) Polyakov loop operator L becomes a 
complex-valued function of the two angles (p 3 and (ps : 


L{p3,ps) ■■= 


1 


exp 


tip3- 




1 

3 

1 

'3 1 


■ IP8- 




e ^ 


-t 26*2^'^ 


COS 


(f) 



C. (III.13) 


C. Our standpoint 

The standpoint of our approach presented in this pa¬ 
per, the first approximation and its improvements, is 
completely different from the other work based on the 
systematic loop calculations in the perturbation theory 
1^ . Although the standpoint of our approach has been 
explained in the previous work Q and this paper is also 
written in this setting, we repeat it below for avoiding 
the misunderstanding on our approach. 

In this paper we aim at a purely non-perturbative ap¬ 
proach in which we look for the initial approximation 
which captures the essential features of the problem in 
question as much as possible at the initial stage, which is 
the spirit of the first approximation. We do not intend 
to do the one-loop calculation in the perturbation theory 
and do not intend to do the systematic loop calculations 
of higher orders, either. Our approach is different from 
[ 5 ^ conceptually in this aspect. We use the terminology 
“one-loop type” to distinguish it from the one-loop in the 
perturbation theory. 

In the first approximation, we take into account only 
the quadratic terms in the fields to obtain the effective 
action (except for the restricted field), which leads to the 
“one-loop type” calculations. It is well known that the 
effective action F obtained from the classical action S 
by the Legendre transform of the generating functional 
of the connected Green functions is equal to the classi¬ 
cal action S plus the additional part represented by the 
logarithmic determinant resulting from the Gaussian in¬ 
tegrations over the quadratic parts. Therefore, the action 
5efi to be calculated by integrating out all fields in the 
first approximation in our paper agrees with the effective 
action F, up to the special treatment of the restricted 
field as explained below. 

The reason of the special treatment of the restricted 
field is as follows. In our formulation, the Polyakov loop 
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operator L[£^] is completely written in terms of the re¬ 
stricted field i.e., L[^] = LlY], If we take the gauge 
(IIL4) or (IIL6), then ‘rf is equivalent to Y as indicated 
in (IIL5) or (IIL7). By integrating out all the fields up 
to the quadratic parts other than the restricted field y, 
i.e., we obtain the effective theory written in terms of 
the the restricted field y alone. This is along the spirit 
of the first approximation mentioned in the above. Then, 
we estimate the Polyakov loop average by the minimum 
of the effective potential obtained from the effective the¬ 
ory. 

The resulting effective theory is identified with the low- 
energy effective theory in the following sense. We use the 
results obtained in the first approximation as the input 
for performing the FRG approach to improve the first re¬ 
sult. In this sense, the first approximation is regarded as 
the initial condition corresponding to the large flow pa¬ 
rameter K at which the FRG analysis start. Or the first 
approximation can be regarded as a preliminary Ansatz 
for solving the flow equation of FRG. The effective ac¬ 
tion with the flow parameter k in the FRG approach 
means that the high-energy modes above k, i.e., p > k 
are already integrated out to obtain the low-energy ef¬ 
fective theory which is valid below k, i.e., p < k. For 
large k at the initial step, the high-energy mode for the 
restricted field is not integrated yet or not to be inte¬ 
grated out, since we identify the restricted field with the 
low-energy modes. This is the reason why the restricted 
field is avoided to be integrated in the first or initial ap¬ 
proximation. While the remaining fields are regarded as 
the high energy modes to be integrated even for large k, 
since they behave as the massive modes and decouple in 
the low-energy region. 

Of course, the above setting is just an approximation 
and cannot be the rigorous treatment and hence this first 
approximation must and will be improved afterwards by 
a systematic method. In fact, we intend to improve the 
first result by the non-perturbative FRG at once (not by 
the systematic order by order loop expansion). This is 
the standpoint of our approach adopted in this paper. 


D. Effective action and effective potential 

In order to obtain the effective potential Ves written 
in terms of the restricted field (similarly for the the 
Polyakov loop L), we perform the functional integration 
over the field variables other than the restricted field 
or the remaining field (massive gluon modes), 
the Nakanishi-Lautrup field ^ (massless scalar mode), 
and the Faddeev-Popov ghost and antighost fields rj, 77 . 

giSetf ['^] 

= J (111.14) 


Then we obtain the effective action SeS- 
SeS 

=£^TT\n[-Dl[G] + M^] + £^Tr\n[-Dl[G] + M^] 

+ IttH-DIIG]] + ^TrH-Dl[G]] 

- TT\n[-Dl[G]] - Tr\n[-Dl[G]] (III.15) 

where is the restricted field variable replaced from 
in a new basis (see Appendix lATl and 

D^[G] := 5 - igG^,, .D^[G] := d + igG^. (III.16) 

Here the first term comes from the integration over the 
remaining field 3^^^, the second term from the integra¬ 
tion over the Nakanishi-Lautrup field ,yV ^ and third term 
from the integration over the FP ghost and antighost 
fields. See Appendix lAl for the details of the calculations. 
Finally, we obtain 

•Seff 

=^^TT\n[-Dl[G] + M^] + ^l^Tr\n[-Dl[G] + M^] 

_ lTrln[-G2[G]] - ^Tr\n[- dI[G]]. (III.17) 

Here we have taken into account only the terms quadratic 
in the fields, in addition to the previous approximation 
where the quantum fluctuation parts y^ are neglected. 
We call the approximations taken up to this stage the 
first approximation. The corrections to the first ap¬ 
proximation are incorporated by using the FRG after¬ 
wards where the explicit dependence on the gauge cou¬ 
pling constant enters in the results. 

In a new basis the restricted field variable G^{x) is 
separated into the background field G^ (cc) and the quan¬ 
tum fluctuation part G^{x), i.e., G^{x) = G^{x) -I- 
Gfj,{x). This separation corresponds to (IHI.9I1 and 
(IHI.IOI) . Therefore, we have the specific uniform (i.e., 
x-independent) background G_^{x) = G^S^q for the re¬ 
stricted field variable in a new basis: 

G^(:r) = Go5;.o + G;.(x). (HI. 18) 

Then, the covariant Laplacian —£)^[G] or —Z)^[G] is 
given 

- {dp T igGp{x)f 

= - d1± 2igGp{x)dp + g‘^Gp{x f ± igdpGp{x) 

= — — 9q ± 2igGodo + g^G^ + [G(a;)-dependent terms] 

= — dg + {ido ± gGo)^ -I- [G(a;)-dependent terms]. 

(HI.19) 

Therefore, the covariant Laplacian has the momentum 
representation at finite temperature: 

— {dp =F igGp)'^ -f {wn ± gGo)'^ + [G-dependent terms]. 

( 111 . 20 ) 
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Note that Go is given by 


SU{2): Go =g"^Ta(^, a = ±1, 

5G(3): Go = ± ^ := (^> 3 ,‘^s), (ni.21) 


where are the positive root vectors of SU{3): 






(111.22) 


IV. SU{2) YANG-MILLS THEORY 


The function F 



A. Existence of SU{2) confinement/deconfinement 
transition 

Symmetries of the SU (2) Polyakov loop average L are 
as follows: 

i) periodicity of 47r in {p: 

L{(p) = L{(p + 4 :Tt) , (IV.l) 

ii) reflection invariance: 


FIG. 1: The plot of F^{p) = f ln[l + _ 

2g-\/h^+j^ cos((^)] as a function of the angle (p for various 
values of M M/T > 0 a,t D = 4. 

The effective action Ses reduces to the effective poten¬ 
tial Ves(p>) written in terms of the background part tp, 
i.e., 


L{p,) = L(-^). (IV.2) 

The Polyakov loop operator is written in a simple form: 

L = cos| e (-1,1], (^e[0,27r). (IV.3) 

According to (111.361) . thus, the Polyakov loop average 
can be calculated at the absolute minimum of the ef¬ 
fective potential. The vanishing Polyakov loop average 
L = 0, i.e., confinement is realized if the effective po¬ 
tential V(ip) has the minimum at (/? = tt. Whereas non¬ 
vanishing Polyakov loop average L 0, i.e., deconfine¬ 
ment is realized if the effective potential V{ip) has the 
minimum at p) ^ t:. 

L = 0 <1=^ ip = TT (confinement: Z{2) symmetric), 

L 7 ^ 0 <1=^ P> ^ (deconfinement: Z(2) breaking). 

(IV.4) 

Symmetries of the SU{2) effective potential V{ip) are 
as follows: 

i) periodicity of 27r in ip: 

Ves{p>) = Ves{p^ + 2 tt), (IV.5) 

ii) reflection invariance: 

VM = V,s{-p>). (IV.6) 

These are the result of Z{2) center symmetry: 


S-eff = V,s{‘p)T-^ J d^-^x (IV.8) 


by neglecting the quantum fluctuation part G, since the 
background part is x-independent. In this approxima¬ 
tion, thus, the effective potential has the momentum rep¬ 
resentation: 


14ff(<p) 

D-l 






--T y] 

2 ^ 

nez,; 

(uJn := 2'KTn). 


nGZ,± 


(27r) 


y\n[{u:n±Tp,f+p^ + M^] 


D-l 


2^ ^ I 

nez.i"' ^ ' 


(IV.9) 


The Matsubara sum (summing over the Matsubara fre¬ 
quencies ojn = 27rTn) can be done to obtain the closed 
form: 


'y ] ln[(a;n -l- T6)^ -f A4^] 

= ln[l -H cosdj, (IV. 10) 

where we have defined 

Ep := ^p^ + M^. (IV.ll) 

Then we have 


VesiP’) = VeB{2TT - p^), (IV.7) 

In fact, these properties are satisfied by the explicit form 
of the effective potential (IIV.9I) obtained below. 


^E / ^^Hi^n + Tp>r+P^ + M^]=T^Fyp>). 

nel,d ) 

(IV.12) 
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Here we have introduced the dimensionless variables p 
and M normalized by the temperature T as 

p ■= p/T, M ■= M/T, (IV.13) 

to define the dimensionless function of the angle 

ip with a parameter M by 

FM =J 

f^{f,ip) :=ln[l + e-^V^^^^-2e-V^^^cos(/3]. 

(IV.14) 


In Fig. [U we have given the plot of F^{ip) as a function 
of the angle ip for various values of M > 0 at U = 4. We 
see that Fj^{ip) is exponentially suppressed Fj^{ip) <C 1 
for large M and vanishes in the limit M —>■ oo uniformly 
in ip, i.e., irrespective of the value of if. Here we have 
used 


pOO 

PMi‘P)=C^D dp p^-'^f^ip^,ip), 
Jo 


Cd :=- 


2D-2^^r{^) 


(IV.15) 


Thus, we obtain the dimensionless effective potential of 
the Polyakov loop average normalized by the temperature 
as 


Vesif) := V,s{f)/T^ = {D- l)F^{ip) - F^{f). 

(IV.16) 

In Fig. [21 we have given the plot of the effective potential 
of the SU{2) Polyakov loop as a function of the 
angle ip for various values ofM :=M/T>0atII = 4. 

At sufficiently high temperature, M = M/T <C 1, the 
gluonic mass M is negligible and the effective potential 
given by 

^{D- l)FoM - F^{ip) = {D- 2)F^{f). 

(IV.I7) 


The high-temperature effective potential has 

the Z{2) breaking minima at = 0,27r, and the ex¬ 
tremum at (/? = TT is a maximum. For D = A, this effec¬ 
tive Mtential reproduces the well-known Weiss potential 


Vw{f) =T^ 

(mod 27r). 


21 




^ I ^ 


(IV.18) 


See the M/T = 0 part of Fig. |2]or Fig. |31 
At sufficiently low temperature, M = M/T ^ 1, 
on the other hand, Fjfj{ip) is exponentially suppressed 
Fj^{if ) <C 1 and the effective potential reduces to 


^ -FoiT)- (IV. 19) 


SUG2) PolyaJcov loop potential 



FIG. 2: The D = 4 effective potential V of the SU{2) 
Polyakov loop as a function of the angle ip for various val¬ 
ues of M := M/T > 0. 


The low-temperature effective potential V/g'^(ip) has the 
Z(/2) symmetric minimum at ip = tt. In fact, the effective 
potential in the sufficiently low temperature is reversed 
to the Weiss potential at sufficiently high temperature 
(See the region M/T > 3 of Fig. |2]or Fig. [3]): 

^ -{D - 2)-^V^^^\ip). (IV.20) 

This indicates the existence of the phase transition 
from the high-temperature deconfined phase to the low- 
temperature confined phase. 

The physical interpretation of this phenomenon is as 
follows. At sufficiently high temperature T ^ the 
gluonic mass M is negligible and hence gluons and ghosts 
contribute equally to the effective potential. Both gluons 
and ghosts are equally responsible for the dynamics in 
the high-temperature deconfined phase. At sufficiently 
low temperature T <C Td, the massive gluons with the 
mass M are not excited by the thermal fluctuation of 
order T < M and do not contribute to the effective po¬ 
tential. The massless scalar mode, ghosts and antighosts 
give the dominant contribution to the effective poten¬ 
tial at low-temperature confined phase. In this sense, we 
can say that the confinement mechanism at finite tem¬ 
perature is the scalar-ghost dominance or unphysical 
mode dominance. 

This result is reliable in the region where the back¬ 
ground part of the restricted field for calculating the 
Polyakov loop average gives the dominant contribution 
to the phase transition compared with the quantum fluc¬ 
tuation part of the restricted field. Therefore, it is appli¬ 
cable to the vicinity of the critical temperature Td- The 
above estimate on the critical temperature is meaningful 
as far as the gluonic mass M does not discontinuously 
change across the transition point, even if it changes as 
the temperature varies. The absence of the discontinu¬ 
ous change of M, namely, continuous change of M across 
the transition point will be shown in an analytical way 
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SU(2) Polyakov loop potential 




A co«fficiait .1’ of the SUf2) effective potential 

Al7 



A coefficient -U of the SU(2) effective potential 

J4 4! 



FIG. 4: The coefficients A 2 and of the SU{2) 

Polyakov loop effective potential Vo{(fi;M) as a function of 
M := M/T at D = 4. 


with the coefficients: 


FIG. 3: The D = 4 effective potential V of 

the SU{2) Polyakov loop for M := M/T = 
0.0,1.0, 2.0, 2.5, 2.6, 2.7, 2.8, 2.9, 3.0, 3.1, (Left) as a func¬ 
tion of the angle ip/n G [0,2), (Right) as a function of the 
Polyakov loop average L = cos ^ G (—1, Ij. 


as well as in the numerical way in the subsequent paper 


«m(p) =21n(l + e a/p"+“") > Q, 
^m(p) = -< 0. 


Cm(p) =■ 


(1 _|_ g-y/p^+M'^-^2 

-y/p'^+M^ _ /^g^-y/p'^+M^ _j_ ^- 2 y/p^+lvP ^ 


12{1 + 


(IV.22) 


Then the effective potential is expanded into the power 
series in the angle variable around (p = tt: 


B. SU{2) Transition temperature and order of the 
transition 


Voip; M) :=Ves,oiT)/T^ = (D - l)F^ip) - Fo/p) 


A. 


A I 2 ,M / x 2 I / ,4 


M\ 

A - 

x2 , ^4,M , 


0,M ' 21 

-b 0{{p - 7rf), 


(IV.23) 


where the coefficients are explicitly given as 


We proceed to estimate the critical temperature and 
determine the order of the phase transition. This is 
achieved by the detailed study of the function Fj^{p). 
The integrand fjQ-{p^,p) of Fj^{ip) is expanded into the 
power series in the angle variable p about p = n (at 
which L = 0): 


/M(P^ t ) =aMip) + bjiiip/ip - 7r)2 -b Cjfy{p){p - 

+ Oiip-7rf), (IV.21) 


^0 M =C^D / dp p^ ^[{D - l)ajfy{p) - ao(p)] 

1 

^^2.M =Cd J dp p°-^[{D - l)bpj(p) - &o(p)], 

1 1 "°° 

^^ 4 ,M =^D J dp p°-^[{D - l)Cj^{p) - Co{p)]. 

(IV.24) 

In the limit M —>■ 0, especially, the coefficient ^ 2,0 is 
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given by 

^^2,0 =Cd{D - 2) y dp p^~'^bo{p) 

poo 

= - Cd{D - 2) / dp p 


°° p-p 

D-2 ® 


(1 + e-P)2 


< 0 , 

(IV.25) 


which recovers the previous result for D = 4 with C 4 = 

2 ^ 


1. W 2 e-P 1 ^ 

—A 2 _n = -^ / dp p = -- < 0. 

Jo (1 + e py 6 


2!'i2.o 


(IV.26) 


Similarly, in the limit M —>■ 0, the coefficient A 4 0 is given 

by 

1 

—214,0 =Cd{D - 2 ) J dp p°~^coip) 

„45_2e-P(l-4e-P + e-2p) 


=Cd(.D - 2) / dp p^ 


12(l + e-P)4 

(IV.27) 


which recovers the previous result for Z1 = 4 [^: 


77^4 ,0 — 

4! TT^ 


‘ .oe-Pil-4e-P+ e-^P) 
dp p^ 


1 


12(1 + e-P)4 


> 0 . 


127r2 

(IV.28) 


Fig.|4]is the plot of the coefficients A^ mI2\ and A^ m/ 4! 
of the SU{2) Polyakov-loop effective potential Vo{(p]M) 
as a function ofM := M/T at D = 4. We observe that 
^2 M < 0 for M € [0, 2.9] and ^2 m > 0 for M > 2.9, 
while ^ > 0 for M G [0, 3.6]. 

There/ore, the phase transition from deconfinement to 
confinement occurs at the temperature Td at which the 
coefficient A^ m changes its signature from negative to 
positive, namely, becomes zero: 


^2,M “ b 

pOO 

^ / dp p^-^iD - l)b^{p) - bo{p)] = 0. (1V.29) 

do 

This condition determines the critical value for the ratio 
Me := M{Td)/Td between the gluonic mass M(Td) and 
the transition temperature Td- For D = 4, we find the 
critical value: 

= 2.9724 , = 0.33643, (1V.30) 

Td M{Td) ^ ’ 

where M may depend on temperature. For instance. 


M(Td) = l.OGeV o Td = 340MeV. (1V.31) 


See Introduction for more information. 


A coefneient 5; of the SU(2) effective potential 



A coefneient of the 5Ul2) effective potential 



Polyakov loop effective potential Vq{L\M) as a function of 
M := M/T at T = 4. 


The above results are rephrased in terms of the 
Polyakov loop average L directly. The angle ip and the 
SU(2) Polyakov loop operator in the first approximation 
is related as 

L = cos ^ —>■ cos if = 2LJ — 1. (IV.32) 

Then the effective potential is rewritten in terms of the 
Polyakov loop average L explicitly. In fact, we have the 
expression for p) in terms of L: 

:= ln[l + e-2Vp=+M= _ {2L^ - 1)]. 

(IV.33) 


In the similar way to the above, the integrand is expanded 
into the power series in L about L — 0: 

/M(P^ t) =dM{p) + hj^{p)L^ + Cj^{p)L^ + 0{L^), 

(IV.34) 

with the coefficients: 

d^{p) =21n(l + e-^^^)>0, 

4g— 

^m(p) =- 


-M 


(P) =- 


_|_ g-\/p2+M2^2 

gg-2-\/ p2+M2 

(1 -p p-\/p'^+M^y 


(IV.35) 


The effective potential is expanded into a power series in 
L (around T = 0): 


VoiL;M) =B, 


0,M ' 


^2,M j^2 _|_ '®4,M 


2 ! 


4! 


+ 0{L^), 


(IV.36) 
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where the coefficients are given by 

pOO 

Bq m =^d / dp p^-^\{D - l)a^(p) - ho(p)] 

1 /*°° . ^ 

^^2,M =C'^ J dp p^-^[{D - l)bj^{p) - 6o(p)], 

1 f°° 

^^4.M =^D dp p^~^[{D - l)Cj^{p) - Coip)]. 

(IV.37) 

Here the fact that only the even powers of L appear in 
the effective potential reflects the center Z{2) symmetry: 
Vo{zL-,M) = Vo{L;M) where z satisfies = 1. Fig.[5]is 
the plot of the coefficients B 2 ^ and ^ of the SU{2) 

Polyakov-loop effective potential Vo{L;M) as a function 
of M ■= M/T at D = 4. 

The coefficient B^ ^ has the same behavior as the co¬ 
efficient ^ 2 M- Both coefficients are negative for 0 < 

M < 2.9 and positive for M > 2.9. The coefficient 

is positive for 0 < M < 3.6. However, the coefficient 
B^ is negative for 0 < M < 1.8, although it is positive 

for M > 1.8. The angle variable can be used for any 
temperature above the critical temperature T^. There¬ 
fore, it is better to use the angle variable p than the 
Polyakov loop average itself L at least in this approxima¬ 
tion. Fortunately, however, the critical value Me = 2.9 of 
^2 M ~ M contained in the positive region which 
is common to both and B^ This allows one to 

study the vicinity of the critical temperature using the 
effective potential Vo{L;M) in terms of L even in this 
approximation. 


C. SU{2) Functional renormalization group 


The result of “one-loop type” for the effective poten¬ 
tial given in the above can be improved using the FRG. 
The flow equation called the Wetterich equation [l^ 
is given for the so-called effective average action F^ 
which depends on the renormalization-group (RG) 
scale k: 


=-STr 


t 







where t is the RG time t := ln-|-, dt '■= -§1 = k-^ for 
some reference scale (UV cutoff) A and R<s>^k is the reg¬ 
ulator function for the field $. Here STr denotes the 
super-trace introduced to include both commutative field 
(gluon) and anticommutative field (quark, ghost). The 
physical result, i.e., the true effective action is obtained 
finally in the limit fc 0, after integrating over the whole 
range of momentum, which is in accord with the origi¬ 
nal idea of Wilsonian renormalization group. See [T^ for 
reviews of the functional renormalization group. 


Although the choice of the infrared cutoff function Rk 
is not unique, it is somewhat similar to the mass term 
with the mass proportional to A:. In the case of massless 
fields, the control of the flow in solving the differential 
equation with respect to k is rather subtle, since the mass 
term originated from the infrared cutoff function Rk dis¬ 
appears eventually in the limit fc 0. This is indeed the 
case of the ordinary Yang-Mills theory, i.e., gluodynamics 
where gluons and ghosts are both massless. 

The Wetterich equation tells us that the “one-loop 
type” expression is the first approximation to the solu¬ 
tion of the flow equation, if the infrared cutoff function 
Rk depending on the flow parameter k is included in the 
relevant “one-loop type” expression. If we choose the spe¬ 
cific gauge in which the color field is uniform, i.e., (1111.41) 
or (|HI.6I) . then the fields relevant in the flow equation 
are the restricted field the remaining field 

the Nakanishi-Lautrup field jV, and the FP ghost and 
antighost r]{x), f]{x), i.e., ‘pf = jV, p, fj). In this 

section, we use the Euclidean formulation. Within the 
approximation adopted in this paper, the flow equation 
is written in the form: 


dtTk =-Tr ■ 




V 


dtRk 


dtVT,M,k, 

(IV.39) 


where VT,M,k is the “one-loop type” part: 

=;^^Tr{ln[G^^ + + S^^Rk]} 

+ ^TT{ln[G^^ + S^^Rk]} 

-TT{ln[G^^ + S^^Rk]}, 

G^^ ■={&p[Y]^P[Y])^^. (IV.40) 

Here, the first, second, and third contributions in VT,M,k 
comes respectively from the integration over the remain¬ 
ing field ^, the Nakanishi-Lautrup field JV^ and the 
ghost-antighost fields 77, p sini. We have used the same 
regulator function Rk for gluons and ghosts up to the dif¬ 
ference due to the tensor structure of gluons. By taking 
the trace over the Lorentz indices, VT,M,k is written as 

VT,M,k :=^^^Tr{ln[G^^ + Rk]} 

- iTr{ln[G^^ -h d^^Rk]}. (IV.41) 

We neglect back-reactions of the temporal component 
Yo on the other spatial components Yj. Assuming an 

expansion around Yj = 0, F^^^ := j^Tk-^ is block- 
diagonal like the regulators. Then the flow equation can 
be decomposed into a sum of two contributions under our 
approximation: 


dtTk = -Tr 


t(2) 


Rk , 


dtRk, 


dtVT,M,k, 

(IV.42) 
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where the gluon regulator Rk,fj,v is the (/i, v) component 
of the block-diagonal matrix in field space: 

Rk^OO Ro^k ^O^kRopt^kiP )5 
Rk^Oj —0 — Rk,jO^ 

Rk,jl =Tjl{p)RT,k = Tjl{p)Zj^kRopt,kT{P^), (IV.43) 

where Tji := Sje — is the transverse projection oper¬ 
ator and i?opt,k(p^) is the (D — l dimensional) optimized 
choice 


Ropt,kiP^) = (fc" - p^)9{k^ - p^). (IV.44) 

The first term in the right-hand side of (IIV.42I) encodes 
the quantum fluctuations of 'fp, while the second one en¬ 
codes those of the other components of the gauge field 
and ghosts. In the present truncation, the second term 
is a total derivative with respect to t, and does not re¬ 
ceive contributions from the first term. Therefore, we 
can evaluate the flow of the second contribution, and use 
its output VT,M,ki"^o) as an input for the remaining flow. 


dtTk 



d^-^p 

(27r)-^~i 


+ dtVT,M,k, 


^(2) 


Rk, 


dtRo,k 

00 

(IV.45) 


In the presence of the infrared cutoff function Rk-, the 
approximate solution of the Wetterich equation is given 
by the “one-loop type”: 


VT,M,k{p) 

D-l. 


rtGZ'^ ^ ' 


+ p"^ + M'^ + Rk{p)] 


1 






n^Z ‘ 


(27r) 


^—\n[{uJn±Tip)^+p'^ + Rk{p)], 

(IV.46) 


where we have introduced the same infrared cutoff func¬ 
tion Rk for all modes. In what follows, we assume that 
the flow parameter k dependence of the gluonic mass M 
is negligible. The VT,M,ki‘p) is obtained by replacing the 
massless remaining gluon X by the massive remaining 
gluon X where the counting of the independent degrees 
of freedom for the massive vector field is different from 
the massless gauge field. The Vr.M.fc is an improvement 
of the expression Vr.fc given for H = 4 as eq.(52) in the 
previous work Q, which is recovered in the limit M —>■ 0. 

Performing the Matsubara sum, we obtain the expres¬ 
sion: 


In the present case, we suppose that the remaining 
gluons exhibit the massive behavior. Therefore, if 
k decreases and becomes smaller than the gluonic mass 
M, the gluonic mass M becomes dominant and persists 
in the limit fc j, 0, even if M can depend on k. In other 
words, the existence of the gluonic mass M guarantees 
a stable renormalization group flow under good control. 
Yet, the ghost is still massless and the control of the in¬ 
frared cutoff effect must be taken into account with great 
case. In this way, the simple “one-loop type” calculation 
based on the massive gluon picture can give rather bet¬ 
ter results than those expected from the naive one-loop 
calculations in perturbation theory. However, this is to¬ 
tally different from the perturbation theory. Indeed, the 
gluonic mass must be generated in the dynamical way, 
which is obviously a nonperturbative result. Thus, the 
justification of the one-loop calculation as a good approx¬ 
imation in the present case comes from the observation 
that the one-loop form in the presence of the infra cutoff 
function Rk is a first approximation to the solution of 
the exact Wetterich equation, provided the gluonic mass 
exists. 


In our treatment, the difference between the three-dimensional 
RG scale kj- and the four-dimensional one k is neglected by 
equating two scales kj- = k just for simplifying the analysis, 
since it is enough for obtaining a qualitative understanding for 
the transition. This is not be the case for obtaining quantitative 
results, see Appendix C of [23| for the precise treatment on this 
issue. 


VT,M,k{p) 

- f -YYjL]r,h , „-2f^^/p^+M^+RUp) 

2 7 (2^)^-i ^ ^ 

- 2e“'^Vp"+M2-kfl;,(p) 

_ iy / A - E- in[l -f g-2/3v'pWRHp) 

- cos ifi]. (IV.47) 


In order to perform the momentum integration, we must 
specify the infrared cutoff function Rk- By choosing the 
infrared cutoff function of the optimal type, Rk{p) = 

— p^)9(kE — p^), we obtain 

VT,M,k{<p) 

= ^^CdT f dp ln[I + 

2 Jo 

_ 2e-/3v'fe^TM^cos(^] 

-b [ dp ln[I -b 

2 Jk 

— cos(/?] 

- IcdT [ dp p^-^ ln[I -b _ 26“^'= cos p] 

2 Jo 

- -CdT / dp ln[I + - 2e-^Pcosp\. 

2 Jk 

(IV.48) 
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Therefore, the k dependent part is separated as 
VT,k{f, M) 

= M) + ^ ^ ^ CdT J dp p°~‘^ 

X {ln[l + e-2/3Vra7^ _ 2g-/3x/F^ 

- ln[l + e-2/3\/?^ - cos p]} 

- \^dT [ dp p^-2{ln[l + - 2e-^'^ cos p] 

2 Jo 

-\ii[l + e-‘^^P-2e-<^Pcosp]}. (IV.49) 

By introducing the dimensionless RG scale k and the di¬ 
mensionless momentum p normalized by the temperature 
T: 

k:=klT, p:=p/T, (IV.50) 

this is rewritten for the dimensionless reduced potential 
defined by 

Vi,{p;M) := VtA^;M)IT^, (IV.51) 


as 


Vf^AM) =VoAM) 


D-1 


-Cd / dpp 
lo 

- 2e" 


D-2 


X {ln[l e-2Vfc^+M^ _ 2e-'^^"+^\os<p] 
- ln[l g-2_ 2e-Vp^+^^ cos p]} 


— t:Cd / dpp "^{ln[l-1-e — 2e '‘cos;/?] 
2 Jo 

(IV.52) 


— ln[l-I-e — 2e ^ cos;/?]}. 


In what follows, we estimate contribution of the k- 
dependent terms to the deconfinement/confinement tran¬ 
sition. We expand the fc-dependent part of AA; M) into 
the power series in :=</? — tt around (/? = tt, using 
cos ip = cos(7r -I- (^) = -1 -I- + 0{A)' 

A ■ A ■ 

tV / iCf\ A I 2i,k,M —2 I 4:,k,Ai ~4 I "Sx 

- —ip H- - —(/? +0((p), 


4! 


(IV.53) 


where the coefficients ^ ^ are given by 
4 


=\,M + Cd I dpp 


D-2 


{D - 1) |ln(l -F _ in(l + g-Vp^+M^) 


{ln(l-He-^)-ln(l-he-P)} 


(IV.54) 


A. 


2,k,M 


=^2,M + 2!C'r> / dpp 


D-2 


D 


_ 1 f g-Vp^+M^ 


,-^/WAA 


2 I (X _|_ g-\/^2+M2)2 (X _|_ g-\/fe2-l-M2^2 


[ e-P 

e-^ ]■ 

y (1 + e-P)2 

(l + e-fc)2 j_ 


(IV.55) 


A 


4,k,M 


D-2 


=^4m + 4!Cd / dpp 
’ Jo 

D-1 j g-Vp=-tM2ji _ 4g-\/p=' + M^ + g-2Vp=+M2] 
2 I 12J1g-\/^7^S7]4 

- (P -5> fc)| 

1 fe-P]! -4e-P-f e-2p] 


12[1 -h e-P]4 


-ip-^k) 


(IV.56) 


Here {p —>■ k) means that p in the preceding term is 
replaced hy k. In the limit M —>• 0, these results reduce 
to those in Appendix B in the previous paper Q. 

Suppose that the reduced effective potential is 

decomposed into two pieces: 

Ktk ■= = A + AV/, (IV.57) 


where the first part A is the (/c-dependent) “perturba¬ 
tive part” obtained essentially by the “one-loop type” 
calculation with the infrared regulator function being 
included, while the second part AV/ represents the “non- 
perturbative part” which is initially zero AV/|fe=A = 0 
and is generated in the evolution of the renormalization 
group. The non-perturbative part AV/ is obtained only 
by solving the flow equation in a numerical way and its 
analytical form is not available at present. 

We expand V^(ip' M) in the power series vci ip = ip — tt: 


AAM) =A, 


A. 


0,k,M 


2,k,M ~2 


A 


4,k,M ~4 


0(ip% 

(IV.58) 


where coefficients A^^ ^ yy can be drawn as functions of 
k and M. Suppose that AFj,((/?; M) is also of the form: 


AT/ ( lC/f\ I ~-2 I ^4,k,lil ~-4 . -Gt 

AV^{ip;M) = aQj^yy + ^^ip + + 0{ip ). 


Then the effective potential has the expansion: 


(IV.59) 


C ~ C ■ 

xyglu©/ aOtx I 2 ,k,M —2 I 4.,k,AI —4 1 ~6\ 

Kff.fe (v?; =<^o.fc.M + +Oip ), 


4! 


(IV.60) 
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where + a„ ^ 

After integrating over the fields other than the re¬ 
stricted field we are lead to the effective action of 


rfc[ro] = i J d^-^x!^^Zo,kd,%{x)d,%{x) + v^^'^^,[%]^ 

[%\ = VT,k m + AFfe [%] ■ (IV.61) 

Then the flow equation is reformulated for AVfc with the 
external input VtX- 


dt{^Vu[%]) = i 


where 


r ( 1 

(2^)3 


dfR 




00 


(IV.62) 


rf [ro] = T-i {Zo.fcp" + d\Vu[%]] . (IV.63) 


Using the specific infrared cutoff function: i?o,fe = 
Zo,k{k'^ — p^)0{k‘^ — p^), which yields 

dtRo,k = [dtZo^kik'^ - P^) + 2^o,fefc^] 6'(/c^ - p^), 

(IV.64) 


we can perform the momentum integration analytically: 


/3dtiAVkl%]) 

2 1 

^3 (2^)2 z-^g^p^d}{VT,k[%] + AUfe[l^o]) + k^ ’ 

(IV.65) 


where ak is the running gauge coupling constant defined 

by 


al 


7-1 9 


(IV.66) 


Z ry — l Z ifn '7'—J 

9k ■- ^o.fc5 I afc ■- 
and rjk is the anomalous dimension defined by 

r]k ■= dt In Zo,fe = -dt In ak- (IV.67) 


The flow equation is simplified for the dimensionless ef¬ 
fective potential and the dimensionless RG scale k as 


dk^Vkin] 


_ (1 -H _ 

^l + ^a2(U^_fe[ro] + AU^[f^o])’ 


(IV.68) 

where all scales are measured in units of temperature. 
The input for the flow equation is just a running gauge 
coupling constant ak, apart from the gluonic mass M 
which is assumed to be independent of k. 

It is shown [3| that the flow equation for the effective 


potential 14 ^“ is reduced to a set of coupled flow equa¬ 
tions for the coefficients in the effective potential (IIV.57I) 


combined with (IIV.58I) and (IIV.59I) : 


1 + ^Vk 


67r2 

1 i^2,k,M + 

(IV.69a) 

1 + Ivk 

ATTak{A^j.j^ + 

67r2 

[1 + k,M + ®2,fc,M)] 

(IV.69b) 

1 + ^Vk 

6[47rafc(A4 fc M -f a4.fc,MF 


[1 + fe2 i^^2, k,M + ®2,fc,M)] 


(IV.69c) 


These equations are coupled first-order ordinary but 
nonlinear differential equations for coefficients ^ jg, de¬ 
rived in Appendix C of Q . We see that this form (IIV.59I) 
is justified as a solution of the flow equation. In fact, it 
is easy to see that k m ~ ^ ^'^d k m ~ ^ 
guaranteed from the flow equation, if the effective poten¬ 
tial has no odd power terms at arbitrary k. Therefore, 
if an initial condition, — Q — a^^j^&tk = K 

is imposed, then = 0 and = 0 are main¬ 

tained for any k G [0, A] by solving the flow equation. 
For the given running gauge coupling constant ak and 
the gluonic mass M, these equations can be in principle 
solved, since the coefficients ^ are given explicitly 
in (IIV.SSP and (IIV.56I) etc.. In practice, however, one 
must truncate the infinite series of differential equations 
(IIV.69cl) up to some finite order to obtain manageable set 
of equations, even if we perform numerical calculations. 

We can understand qualitatively why a 2nd order 
phase transition from the deconfinement phase to the 
confinement phase can occur by lowering the tempera¬ 
ture. At a certain temperature T, the flow starts from 
the “one-loop type” result without non-perturbative part 
at fc = A I. Therefore, a^i,j^ = Oa.tk = A:s>l and 
hence — ^2 k m atfc = A^ 1. Moreover, we 

assume + ®4,fe,M > 0 for 0 < A: < A, 

as a necessary condition for realizing a 2nd order tran¬ 
sition. Otherwise, we must consider the higher-order 
terms, e.g. This assumption allows us to ana¬ 

lyze just one differential equation for obtaining qualita¬ 
tive understanding. Then the right-hand side of (IIV.69bl) 
is negative for any k G [0,A], since the running cou¬ 
pling constant ak is positive and 1 -I- iryfc is positive, see 
Fig.6 and Fig.7 of Consequently, 02 ,fc started at zero 

becomes positive 02 ,fe > 0 just below A and increases 
(monotonically) as k decreases. See Fig.8 of 0 for the 
massless case M = 0. Thus, the flow always moves in the 
direction enhancing confinement. If we consider the low 


This assumption is assured to be true by numerical calculations 
of the full effective potential [2^ . [ 2 ^ , as reproduced in the pre¬ 
vious section. 
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temperature T <Td where ^2 fe m ^ then the final re¬ 
sult is always C'2j,^>0atfc = 0. Even if we start from 
^2 fc M < 0 at temperature T slightly above T^, T > Td 
where = ^ 2 ,k,M ® 2 ,fc,M < 0 at ^ ^ ^ 

it may happen that 0 ' 2 ,fc,M = \,k,M + ^ 2 ,k,M > 0 at 
k = 0. But, the ERG improvement does not change the 
above conclusions in an essential manner. Thus, we con¬ 
clude without the detailed numerical calculations that 
the above Td gives a lower bound on the true critical 
temperature Tcj since the flow evolves towards enhanc¬ 
ing the confinement, under the assumption that M does 
not change so much along the flow. 

Finally, we mention the pressure in the ERG. The pres¬ 
sure Pk{T) at the flow parameter k is defined in the low- 
temperature confined phase by Pk = “ Train = 

= -Co.fe.M = -\k,M - %,k,M- 
tion is Oq jg = 0atA: = A:^lor = ^o,fe,M 

at A: = A ^ 1. Then, the flow of the pressure Pk 
is determined from the behavior of Og ^ ^ governed by 
(IIV.69al) . In the low-temperature confined phase, the 
right-hand side of (|IV.69ap is positive, since k m = 
^2 fc M fc M ^ yields the positivity of the 

derivative, ^ m > 0, that is to say, is mono- 

tonically decreasing in decreasing k. Therefore, Pk is 
monotonically increasing in decreasing k due to the ERG 
improvement and finally reaches the largest value for the 
true pressure P at k = 0. Therefore, the true pressure 
can be positive, even if the initial pressure is negative 
in the initial approximation. See Fig. [51 The Pk{T) in¬ 
creases more rapidly at the temperature T closer to the 
critical temperature Td than that at the low temperature 
T <C Td, since C 2 k,M^'^) ~ ^ 2 ,fc,M(^) ® 2 ,fc.M(^) 

comes smaller as the temperature is closer to the critical 
temperature, leading to the larger value of the derivative 
^k^okM according to pV.69al) . In this way, the ERG 
will improve the positivity violation of the entropy in 
the first approximation near the critical temperature in 
the low-temperature confined phase. See Fig. 0 This 
tendency agrees with the two-loop improvement of the 
one-loop resul t 15811 and and is consistent with the other 
ERG analysis [3l|| . 


D. SU{2) Pressure and entropy 


Our effective potential is given by 

Kff(T) 

={D - 1)TCd f dp ln(l -H 
Jo 

— 1"^ cos ip) 

pOC 

-TCd / dpp^“2ln(I-he-2p/^-2e-P/^cos(^). 

Jo 

(IV.70) 


This should be compared with the other effective po¬ 
tential obtained in One finds that the 


effective potential contains two extra terms: 

In fact, the difference is given by 

={D-1)TCd [ dpp^-^ln(l-e-Vp"+“'/'^) 

Jo 

pOO 

- TCd / dp p^-^ ln(l - 6“^/^), (IV.71) 

Jo 

where the first and second terms come from the neutral 
(or the diagonal) component for the gluon and ghost re¬ 
spectively. (The reason of this difference is explained 
below.) However, both effective potentials give the same 
value for the critical ratio M/Td- (Notice that the dimen¬ 
sionless effective potential Vesip)/T^ is written in terms 
of the dimensionless quantity M/T.) This is understood 
as follows. The critical value Td is determined by the 
Polyakov loop average so that it separates the confined 
phase T = 0 in low temperature T < Td and the decon- 
fined phase L ^ 0 in high temperature T > Td- The 
Polyakov loop average L is calculated as T = cos 
using the value Pmin which gives the minimum of the ef¬ 
fective potential t4gf: the derivative is equal to zero at 

V^min • 

Kff(Tmi„) := = 0. (IV.72) 

But, the extra terms (IIV.71I) do not depend on p, and 
therefore do not change the location of the minimum 
Train- Thus we obtain the same value (^min) which 
implies the same value for the Polyakov loop average 
L = cos for a given ratio T/M. Thus we obtain 
the same critical value of the ratio M/Td for the two 
different effective potentials. 

However, the two effective potentials have the different 
minimum values even at the same value (^min due to the 
:/5-independent extra terms (IIV.7II) . 

The pressure P is defined through the temperature- 
dependent minimum value of the effective potential: 


P(T) := -VesiT = Traia(T)). (IV.73) 

The extra terms give the T/M dependent shift for the 
pressure P. In our case, we have 

~ ~(d^ ~ l)^M(Tmin) + Po(Tmin)j (IV.74) 

where Fm{t) is defined by (IIV.14I) and (I IV. 151) . In the 
low-temperature confined phase, leads to the 

pressure: 

P^ =- 2{D - 1)TCd [ dp ln(l -f 
Jo 

pOO 

+ 2TCd dpp^-'^\n{l + e-P/^). (IV.75) 

Jo 

In the low-temperature limit M := M/T —>■ oo, the mini¬ 
mum of the effective potential is given at Train = t’", which 
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SU(2) Pressure rescaled by T‘^ 

p/7-4 



FIG. 6: The pressure P (rescaled by T'^) as a function of T /M 
in the low-temperature confined phase for SU (2) and D = 4. 
The top line is our pressure P^/T*, the bottom line is the 
RSTW pressure middle line denote the 

difference of the two pressures p^ jT'^ — pRSTWyy4 jjgj-g ^j^g 
critical value is Td/M = 0.33643. 


yields the positive value: 

^ ^ -{D- l)Foo(7r) -h Fo{tt) = Fo(7r) 

^ 2 ^ ^ g-p) ^ ^ 0.191909 {D = 4), 

(IV. 76) 


which is the behavior for the gas of free relativistic mass¬ 
less fermions with 2 internal degrees of freedom. See 
Fig. m In the high-temperature limit M := M/T —>■ 0, 
the minimum of the effective potential is given at ifmin = 
0, which yields the positive value: 


pK 

pD 


-{D- l)Fo(O) + Fo(0) = (2 - D)FoiO) 
' dp 


-2x2 


-1 


27r^ 


ln(l — e 


= -2 X —TT^ = —TT^ ~ 0.438649 (0 = 4), (IV.77) 


45 


45 


which is the behavior for the gas of free relativistic mass¬ 
less bosons with 2{D — 2) internal degrees of freedom. 

The difference between ours and RSTW is pmin inde¬ 
pendent, and depends on the ratio M := M/T alone for 
a given dimension D: 


pRSTW 

pD 


pK 
pD 


D 


=Cd 


[D-l) / dgg"-^ln(l 




dq ln(l - e-?) 


(IV.78) 


For D = 4, the difference is given by 


pK pRSTW 2^ 

r4 “^2 L' 


3 / dqf ln(l - 

Jo 


dq f ln(l - e-?) 


(IV. 79) 


In the low-temperature confined phase T < Td, both 
pressures p^/T'^ and pR-STW^p4 increase monotonically 
in M = M/T or decrease monotonically in T/M. (This is 
not the case in the high-temperature deconfined phase.) 
The difference (I IV. 79(1 is monotonically increasing in 
M = M/T or monotonically decreasing in T/M in both 
phases. In the low temperature M := M/T ^ 1, the 
difference is positive: 

pK pRSTW 1 /-oo 

=—TT^ ~ 0.109662 for D = 4. (IV.80) 
90 ^ ' 


In the high temperature M := M/T <C 1, incidentally, 
the difference is negative: 


pK pRSTW 
2~’4 


1 

27r^ - 

1 


dq f ln(l-e-«) 


= - —TT^ ~ -0.219325 for D = 4. 

45 

(IV.81) 


The RSTW pressure pi^STW violate 

slightly the positivity before reaching the critical tem¬ 
perature as the temperature is increased in the low- 
temperature confinement phase. Indeed, the positivity 
is maximally violated pRSTW^p4 _ _o.00161342 at the 
critical temperature Td/M = 0.33643, while our pres¬ 
sure is positive /T'^ = 0.0217697 even at the critical 
temperature where the difference is 0.0233832. Thus, our 
pressure in the initial approximation remains positive 
in the low-temperature confined phase, in sharp contrast 
to the RSTW pressure at one loop. 

This difference is understood as follows. Our mass 
term does not agree with the mass term introduced in 
[53 even after fixing the gauge as performed in section 
III. The two mass terms have different independent de¬ 
grees of freedom. Our mass term is written in terms of 
the remaining field alone which are charged. This fact 
for the independent degrees of freedom for the remaining 
field yields the ghost fields different from those in [^ . 
The authors of introduce the neutral ghost associated 
with the neutral or diagonal gauge field, in addition to 
the charged ghosts associated to the charged gauge fields. 
Notice that the neutral ghost gives the negative contri¬ 
bution to the pressure, while the charged ghosts give the 
positive contribution to the pressure. However, in the 
gauge fixed version of our formulation, the neutral ghost 
is prohibited to be included and is not introduced even 
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SU(2) Pressure 

?IM* 



FIG. 7: The pressure P (rescaled by M'*) as a function of T/M 
in the low-temperature confined phase for SU{2) and D = 4. 
The top line is our pressure the bottom line is the 

RSTW pressure and the middle line denote the 

difference of the two pressures P^ /M^ — /M'^. Here 

the critical value is Td/M — 0.33643. 


densations. This suggests that I4ff(L) is insufficient to 
examine the minimum value near the critical temperature 
in the low-energy confined phase and is to be replaced 
by the simultaneous effective potential of $ 

and L, since the non-vanishing temperature-dependent 
condensate $ defined by dEH) will lower the vacuum en¬ 
ergy to give a different temperature dependence for the 
pressure. The gauge-invariant gluonic mass M obtained 
from $ could be related to the glueball mass, see [ 23 - 
Taking such dynamical degrees of freedom for glueballs 
into consideration is expected to eliminate the artifact 
of considering L alone to recover the positive entropy, 
i.e. monotonicity of the pressure. Indeed, there exist 
other works suggesting that the glueball degrees of free¬ 
dom reproduce the expected thermodynamic behaviors 
of the Yang-Mills theory, see e.g., [13,0, llll • The result 
of the effective potential will be reported in a 

subsequent work. 


after fixing the gauge for the remaining fields. This is 
because the remaining field .^(a;) is required be orthog¬ 
onal to the color field n(x), namely, the defining equation 
(II) ■ n(x) = 0 in section II must be imposed, and 

hence the ghost fields associated with the remaining fields 
must be also orthogonal to the color field, which means 
that the ghost fields must be charged, in other words, 
the neutral ghost, i.e., the component of the ghost field 
which is parallel to the color field must be vanishing. (See 


111, section 4.6 and Appendix E for S/7(2), and section 


5.8 and Appendix H in for SU(N).) Therefore, the 
remaining fields and the associated ghost fields have less 
independent degrees of freedom than those in [s^ . In 
fact, this difference avoids the violation of positivity of 
the pressure in the low-temperature confined phase. 

In view of these, our initial approximation is relatively 
good. Of course, the first approximation is not enough to 
derive all essential aspects of the finite temperature Yang- 
Mills theory. Some results are to be improved to avoid 
the artifacts of the initial approximation. For instance, 
we consider the entropy density S defined by 


SiT) := . (IV.82) 

We observe that the pressure P is increasing at small 
temperature and hence the entropy S(T) is positive. 
However, as the temperature is increased, the pressure 
changes its monotony and begins to decrease, indicat¬ 
ing that the entropy becomes negative in the region 
Tq < T < Td before reaching the critical temperature 
Td, which was pointed out in See Fig. [71 Here 

To/M = 0.286 for our case and Tq/M = 0.266 for RSTW. 

The positivity of the entropy density S(T) is equal to 
the monotonic increase of the pressure P{T) with respect 
to the temperature T, which means that the minimum 
I4ff(Tmin) of the effective potential I4ff(T) is monotoni- 
cally decreasing and the bottom becomes deeper as the 
temperature increases. The vacuum energy is further 
lowered by existence of more non-vanishing vacuum con- 


V. SU{3) YANG-MILLS THEORY 

A. Existence of SU{3) confinement/deconfinement 
transition 


Symmetries of the SU{3) Polyakov loop operator L are 
as follows: See Fig. jS] 

i) periodicity of 47r in the (/?3 direction and 4-\/37r in 
the ips direction: 

=L{ip3 -b 47r, (fis) = L{ip3, ips + 4-/37r), 
=ReL{ip3 + 47r, tps) 

=ReL{(f3, (fs + 4-/37r), 

=ImL{(p 3 + 47r, (ps) 

=/mT(v?3, (/?8 + dv^Tr), (V.l) 

ii) reflection symmetry: 

L{P3,<P8) = L{-(p3,(ps), 

=>ReL{ip3, (fis) = ReL{-ip3, ips), 

ImL{(p3, (fis) = ImL{-(p3, (ps), (V.2) 


L{P3,V8)* = L{p3,-p8), 

ReL{p3, ps) = ReL{p3, -ps), 

ImL{p 3 , ps) = -ImL{p 3 , -ps), (V.3) 


L{P3,P8) 

ReL{p3,ps) 

ImL{p3,ps) 


iii) global color symmetry: 

L{p3,Ps) = L{-P3,-P8), 

^ReL{p'3, (/jg) = ReL{p3, ps), 

ImL{p'^, Ts) = -ImL{p 3 , ps), (V.4) 
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FIG. 8: 3D plot and contour plot of the SU{3) Polyakov loop 
as a function of the two angles ifiz/n and (Left) Real 

part, ReL, (Right) Imaginary part, ImL. 


where (</? 3 ,<P 8 ) i® obtained from by a rota¬ 

tion of angle ±7r/3: 


¥23 


cos ? ± sin ^ 


¥^3 


’ 1 v^’ 

2 2 


¥23 

_¥28_ 


1 

-H 

0 

0 

Ol 


Ts 


1 

_“r 2 2 . 


¥28 


(V.5) 


ii) charge conjugation invariance: 

-•^s) = Ves[-^P3,'^s), (V.9) 


iii) global color symmetry: 


v;ff(v53, = v;ff((p 3 , (/7g), (v.io) 

where ((/Sg, (/ 3 g) is obtained from ((/? 3 , (/jg) by a rota¬ 
tion of angle ± 7 r/ 3 : 


¥^3 


cos ^ 

zb sin ^ 


¥^3 


’ 1 Vs’ 

2 2 


¥^3 

¥^8 


=Fsin| 

cos ^ 




1 

“T 2 2 


Ts 


(v.ii) 


The transformation (IV.1I|) is equal to 

r_V3 


/ I 2^3 

¥^3 = •^ 


i (^3 - #¥^8 


¥=8 = 


2 ¥>3 + ^¥^8 
+ ■^¥’3 + ^¥’8 


which leads to 



(V.12) 


(V.13) 


The transformation (IV.5I) is equal to 


/ - ; 2 


¥=3 = 


+ #¥28 


i¥23 




which leads to 
-2 , 


¥28 = 


V3 


¥>^ + ^8 = 


-¥2^ + ^8 = 


. ¥28 = 


-¥23 ■ 


-(¥^3 + 75¥’8 


-#^3 

-#¥23 






-¥^3 + ^¥^8 

#3 + ;)|¥’8) 


) 


Vs 


¥28 


5 ¥28 
^¥28 


(V. 6 ) 


(V.7) 


It is easy to see that the Polyakov loop operator (IIII.13I) 
respects all the symmetries i), ii) and iii). 

Symmetries of the SU{3) effective potential 
Veff(¥ 23 , ¥28) are as follows 

i) periodicity of 47r in the ips direction and 47r/-\/3 in 
the ifs direction: 

Vesi^PSi^s) = 14fi(v23 -l-47r,(/j8) = Kff(¥23,¥28 -l-47r/v^), 

(V.8) 


From the above symmetries, it is sufficient to con¬ 
sider the background-field effective potential VgsiPs, Ps) 
in an equilateral triangle OAB of side length 
with the three vertices at O : (0,0), A : (27r, 27r/-\/3) 
and B : (27r, —27r/-\/3) in the {(p 3 ,(ps) plane. The ef¬ 
fective potential Ves{p 3 ,pn) is invariant under rotations 
which leave this equilateral triangle invariant. The ver¬ 
tices O, A, B of the triangle and its center G located at 
(47r/3, 0) are always extrema of the effective potential 
Therefore, we consider the effective poten¬ 
tial at the points 0,A,B^ and G and their vicinity. See 
Fig. and Fig. [TUI 

At high temperature M/T <C 1, the minima of the ef¬ 
fective potential are realized at the points O, A, B (the 
vertices of the triangle) which correspond to the non- 
vanishing Polyakov loop average L = 1, e®3'^ re¬ 

spectively, i.e., deconfinement. Choosing one of them 
spontaneously breaks the center symmetry. The ef¬ 
fective potential takes the maximum at the center point 
G which corresponds to the vanishing Polyakov loop av¬ 
erage L = 0. 

At low temperature M/T ^ 1, the minimum of the 
effective potential is realized at the center point G with 
vanishing Polyakov loop average L = 0, i.e, confinement. 
The points O, A, B (the vertices of the triangle) are the 
maxima which correspond to the non-vanishing Polyakov 
loop average L = 1, respectively. 

The above statements are summarized in the following 
table. 
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SU C3) Pol 3 rakciv bop potentbl at M/T dD 



FIG. 9: Plot of the D = 4 effective potential V of the SU{3) 
Polyakov loop as a function of the two angles (fis /tt and tps /tt 
(L eft) at M := M/T = 0, (Right) at M := M/T — 3.0. 


SUfj)P<^yak<n'looppoteatialat M/T=0 



FIG. 10: Gontour Plot of the 0 = 4 effective potential V of 
the 51/(3) Polyakov loop as a function of the two angles ifia/n 
and (Left) at M := M/T = 0, (Right) at M := M/T = 

3.0. 
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By representing the trace explicitly as the sum over 
the Matsubara frequencies and the integration over the 
spatial momentum, we obtain the effective potential 
I 4 ff ((^ 3 , (/jg) in terms of the two angles 933 and (/Jg: 


The effective potential 14 ff(:p 3 , (/Jg) is explicitly written 
as 


f4ff(<P3,<P8) 
D- 1 


+ In ■ 




nGZ.i 


d^-^p 


ln[(a;„ ± TpsY + + M^] 


±rI -ifs + —ifs 


+ p'^ + m'^ 


In ■ 


± r I -ip3 - —ifs 


+ p'^ + M'^ 




p 


Vee{T3:T8)=^^^Tj2Y. I (27r)^-i 


X ln[(a;„ + ■ pY + p'^ + M'^] 

-YT.T. ' 


P 

2^ ^ ^ J (27r)^-i 


X ln[(a;„ + + p\ (V.14) 






2 ^ J (27r)^-i 

nez.i-^ ^ > 


h\[{uJn±Tip3Y +P^] 


In ■ 


In ■ 


Wn ± r I -ip3 + 


Wn±TI 2 /^ 3 - 


+ P" 


P 


. (V.15) 
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After performing the sum over the Matsubara frequen¬ 
cies uJn = 27rTn, thus, we obtain the effective potential 
for the Polyakov loop average as 




D 


={D-l) 


Fm{^3) + F, 


M 


7^3 • 


73 


-<P8 


+ F-, 


M 


2^3 


73 


-‘Ps. 


Fq{(Pz) + -Fb I 2'^3 + j + -Fb 


T<F3 


73 \ 

— 

(V.16) 


At high temperature, M = M/T <C 1, the mass M is 
neglected and the effective potential reduces to 




^o(V53) + ^0 I ^<^3 + ^‘F8 


+ -Fb I 2‘F3 - —‘/?8 


(V.17) 


For D = A, the effective potential reproduces the well- 
known SU(3) Weiss potential This potential has 

degenerate minima on the vertices of the basic equilateral 
triangle, leading to a deconfined phase with the sponta¬ 
neously broken Z(3) center symmetry. 

At low temperature, M = M/T ^ 1, on the other 
hand, {(p) is exponentially surpressed ((/?) <C 1 and 
the effective potential reduces to 


D 


V,^riT3,T8)/T 


Foips) + ^b I 


Fo(^lps-^Tsj 

(V.18) 


The effective potential at the sufficiently low tempera¬ 
ture is reversed to the Weiss potential at sufficiently high 
temperature: 

Ts) ^ -{D - (^8). (V.19) 

Therefore, the effective potential has the absolute min¬ 
imum at the center G of the triangle leading to a 
center symmetric confining phase. Thus there must ex¬ 
ist a phase transition at a certain critical value of Td/M 
between the high-temperature deconfined phase and the 
low-temperature confined phase. 


B. SU{3) Critical temperature and order of the 
phase transition 

The absolute minimum of 14ff(v^3, Ts) lies on the ps = 
0 axis up to the discrete rotations by the angle tt/S for 


all temperature: 


Ves(p 3 , 0)/T^ ={D - 1) [F^ips) + 2F^ 
- Foip3) + 2Fo{^) 


(f) 


where 

FiiiiTs) 




{2tt) 


D-1 


(V.20) 


/M(P^‘F3) + 2/^(p^^)] . (V.21) 


Fig. [mis the plot of the Polyakov-loop effective potential 
Vo{p 3 , 0; M) := 14ff(<P3, 0)/T^ at (ps = 0 as a function of 
ps for various values of M/T in D = A dimensions. 

The integrand fpfip^jPs) is expanded into the power 
series in ps about ps = Att/S at which L — 0: defining 
a := ps — 47r/3 


Ts) =c® (p) + c7 {p)a + (p)cr2 + (p)cr3 


+ c2Hp)7 + 0(a7, 


(V.22) 


where the coefficients c71 (p) are given in Appendix |B] 

Then the expansion of the integrand fpdip^jPs) + 
into the power series of a is given by 


/M(P^ + 2/J^^(p^ y) 

=h%^ (P) + (P)<^ + (Py^ + 

+ /i2^(p)7 + 0(a7, (V.23) 


where the coefficients h^^\p) are given in Appendix |BJ It 
should be remarked that the linear term in a disappears 
finally = 0. Accordingly, the effective potential 

has the power series expansion in cr: 


7o(p3,0;M) ■.= Ves,o{p3,0)/T^ 

A 


„ , ^2,M_2 , ^3,M_3 


4! 


-7 + 0 {a^), 

(V.24) 


where the coefficients A^ {n = 0,1, 2,...) are given by 

pOO 

^u,M =n'-CD / dp p^-^D - l)h7^(p) - h!F^{p)\. 

J 0 

(V.25) 


In the limit M —)> 0, especially, we find 

poo 

A„,o =n!Cz3(F-2) / dp p^-^hi^\p), (V.26) 

Jo 


which reduces for _D = 4 to 

^2,0 = ~ o ^ ^3,0 — 7 ;~ — 0.15916... > 0, 

o Ztt 

A 4.0 =-^ = 0.22797... > 0. (V.27) 
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SU C3) Polyakov bop potential at ^ HD 



FIG. 11: The D = 4 effective potential V of the SU{3) 
Polyakov loop at ips = 0 as a function of an angle ips/ir € 
[—1, 3) for various values of M := M/T. 


Note that there is no linear term of a in the effective 
potential. Fig. [T^ is the plot of >l 2 m> ^3 M’ ^4 m 
the SU(3) Polyakov loop effective potential as a function 
of M ■= M/T at ZD = 4. 

In order to discuss the order of the phase transition, 
we consider the effective potential (the Landau function) 
of the form: 

V(a) = —ha — -aa^ + -ba^ + -ctr^. (V.28) 
2 3 4 

The extrema are obtained by solving the equation: 

V'{a) = —h — aa -\- ba^ + ccr^ = 0. (V.29) 

We restrict our consideration to a case in which h is neg¬ 
ligible, i.e., h = 0, which is the case for the SU{3) Yang- 
Mills theory. The stationary points of V (a) at h = 0 are 
given at three values of a, i.e., a = 0, cr+, cr_: 

—b + ^6^ -I- 4ac —b — 

a = :=--. 

(V.30) 

At the stationary points a = 0,(T+,cr_, the effective po¬ 
tential has the values: 


V{a = 0) = 0, 
V{a±) = - 


b^ + Qab'^c + =F b{b‘^ + 4ac)^/^ 


24c3 


(V.31) 


The first order transition occurs when the two minima 
give the same value of the effective potential (free en¬ 
ergy), namely, the condition V{a = 0) = V{a±) is satis¬ 
fied: 


2 5(T)2 
9^- 


SU(3) Polyakov loop potential at ^3=0 
VT^ 



v/r’ 



FIG. 12: The D — A effective potential V of the 
SU{3) Polyakov loop at yjg = 0 for M ~ M/T = 
2.65, 2.70, 2.75, 2.76, 2.80, 2.90, (Left) as a function of an an¬ 
gle (fis/ii G [—1,3), (Right) as a function of the Polyakov 
loop average L — ^ [l-|-2cos(^)] G (—1/3,1], normalized 
as F(L = 0) = 0. 


at which the global minimum experiences a discontinuous 
jump. This condition determines the critical temperature 
Td- The first order phase transition is induced by the 
cubic interaction ^ba^. When 6 = 0, the condition (IV.321) 
reduces to a = 0 as long as c > 0. This is nothing but the 
condition for the second order phase transition, which is 
indeed the SU{2) Yang-Mills theory. 

The confinement/deconfinement phase transition in 
the SU{3) Yang-Mills theory described by the effective 
potential (IV.241) is of the first order. See Fig. [121 The first 
order phase transition for confinement/deconfinement in 
the 51/(3) Yang-Mills theory is induced by cubic inter¬ 
action L^, which occurs when the condition is satisfied: 

^ 2 ,M — 3 (V.33) 

See Fig. [14] for the plot of mA^ ^ — 5(^3 ^)^ for the 
SU{3) Polyakov-loop effective potential as a function of 
M := M/T at ZD = 4. This shows that the condition 
(IV.331) for the first order transition is indeed satisfied at 
M = 2.75, which is greater than M = 2.6 at which the 


a(r) = 


(V.32) 
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A coefficient A2 of the SUf3} Polyakov loop potential at =0 
A- 



A coefficient A; of the SUf3) Polyakov loop potential at =0 
•^3 




A coefficient £2 for the SU(3) Polyakov loop potential st.fs=0 



A coefficient £3 for the SU(3) Polyakov loop potential at =0 

S 3 



A coefficient ^4 for the SU(3) Polyakov loop potential at =0 
54 



FIG. 13: The plot of ^ 3 ,M) S'lid for the SU{3) 

Polyakov loop effective potential as a function of M ~ M/T 
at D = 4. 


FIG. 15: The plot of B^ j^, and B^j^, for the SU{3) 

Polyakov loop effective potential as a function of M := M/T 
at D — 4. 


A2A4 -(1/3) A3 for the SU{3) PoK'akov loop potential at >^8=0 

^ 2^4 ^ V '' 3)^3 ^3 



FIG. 14: The plot of 412 ^^44 ^,;^ — |(4l3 for the SU{3) 
Polyakov loop effective potential as a function of M := M/T 
at D — 4. 


second order transition A 2 ^ = 0 would be realized. In 
other words, when the temperature is decreased starting 
from the high-temperature deconfined phase, the phase 
transition to the low-temperature confined phase occurs 
at a temperature lower than the expected temperature at 
which the coefficient changes its signature from negative 
to positive. 

The transition can be observed by the effective po¬ 
tential directly written in terms of the gauge-invariant 
Polyakov loop average. Due to the existence of the center 
symmetry, the real-valued effective potential must have 
the general form: 


Vo{L,0;M) =a, 


O.M ' 


^2,M 


L*L + 


^3,M 


Re{L^) 


+ ^{L*LY + 0[L^), 


(V.34) 


where L is complex-valued in general. In fact, we find 
that the effective potential of this form reflects the ^(3) 
center symmetry: Vo{zL,0] M) = Vo{L,0-, M) where z 
satisfies zz* = I and = 1. 
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At </38 =0, the Polyakov loop operator L is real-valued 
and given by 


L 


1 

3 


1 - 1-2 cos 



(V.35) 


At vjg = 0, the angle is related to the real-valued 
Polyakov loop operator as 


which yields 

costps = 2 cos^ - t = 2 { ^^2 


For tps = 0, therefore, the effective potential has the 
form: 


B 


B. 


-'o.M 

B 


2 ! 


3,M^3 


3! 




(V.38) 


If the first order transition is induced by the cubic term 
3 ',^ Li ^, then the transition from deconfinement to con¬ 
finement occurs at the temperature Td at which the con¬ 
dition is satisfied: 

^2,M — 1^4,m- (V.39) 

This condition determines the value of the ratio be¬ 
tween the transition temperature Td and the gluonic mass 
M(T) which may depend on temperature. See Fig.fTSlfor 
the plot of ^2 Xf, 7^3 and B^ for the SU (3) Polyakov 

loop effective potential as a function of M := M/T at 
D = A. 

We find that the ratio between the transition temper¬ 
ature Td and the gluonic mass M(T) is given for Z? = 4 
by 


M{Td) 

Td 


2.75 


Td 

M{Td) 


0.364. 


(V.40) 


SU(3) Pressure =-Vsii{ip^ = —7:, ^P8=0) rescaled by T‘^ 

VIT* 



T/M 


P/M”* 



T/M 


FIG. 16: The pressure P as a function of T/M in the low- 
temperature confined phase for SU{3) and P = 4. (Left) 
rescaled by T*, (Right) rescaled by AT*. Here the critical 
value is Td/M = 0.364. 


C. SU{3) Pressure and entropy 

The pressure P is defined through the temperature- 
dependent minimum value of the effective potential: 

P(T) := = iprainiT)). (V.43) 


For instance, 

M{Td) =0.8GeV o = 291MeV. 


(V.41) 


In the first approximation, we have 


This should be compared with the zero-temperature re¬ 
sult: 


M{T = 0) = 0.8 - l.OGeV. (V.42) 

This result should be compared with the work : for 
P = 4, the gluonic mass parameter m = 510MeV was 
obtained from fits of SU{3) lattice data for the gluon 
propagator in the Landau gauge at zero temperature, 
which gives the estimate on the transition temperature 
Te = 185MeV. 


P/T^ = = 


D 


= -{D-l) 


F^(^r) + 2F^ 




+ 2Fo 




(V.44) 


where Fm{t) is defined by (IIV.14I1 and (llV.lbp . See the 
first figure in Fig. 1161 In the low-temperature limit M := 
M/T —>■ 00 , the minimum of the effective potential is 
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given at which yields the positive value: 


P/T^^Fo(-7r) + 2Fo -TT 


-^TT^ ~ 0.316802 {D = 4). 
405 ^ ' 


(V.45) 


In the high-temperature limit M := M/T —>■ 0, the mini¬ 
mum of the effective potential is given at = 0, which 
yields the positive value: 

P/T^ -(D-2)3Fo(0) 

^-6 X — TT^ = -^TT^ ~ 1.31595 (£1 = 4). (V.46) 

45 15 

We observe that the pressure P{T) remains positive at 
any temperature in the low-temperature confined phase 
T < Td — 0.364M. However, we find that the simi¬ 
lar phenomenon to the SU{2) case occurs also in the 
SU{3) case. The pressure P{T) is increasing at small 
temperature and hence the entropy S(T) is positive. 
However, as the temperature is increased, the pressure 
changes its monotony and begins to decrease, indicat¬ 
ing that the entropy becomes negative in the region 
Tq < T < Td before reaching the critical temperature 
Td- Here Tq/M = 0.29904 for the SU{3) case. See the 
second figure in Fig. 1161 Therefore, we need the improve¬ 
ment for SU{3) case similar to the SU{2) case, which will 
be reported in a subsequent work. 


VI. CONCLUSION AND DISCUSSION 

In this paper we have shown the existence of the con¬ 
finement /deconfinement phase transition at a finite tem¬ 
perature Td in SU{2) and SU{3) Yang-Mills theories by 
calculating the effective potential of the Polyakov loop 
average. The key ingredient to derive the phase transi¬ 
tion is the introduction of a dynamically generated glu- 
onic mass M in the reformulated Yang-Mills theory which 
allows one to introduce the gauge-invariant mass term for 
a specific gluonic degree of freedom. The transition tem¬ 
perature Td is estimated as the ratio to the gluonic mass 
M. The transition is continuous for SU{2) and discon¬ 
tinuous for SU{3). 

The existence of the gluonic mass has been estab¬ 
lished at zero temperature and has already played the 
very important role in understanding quark confine¬ 
ment at zero temperature. The existence of the glu¬ 
onic mass across the transition temperature enables 
one to easily understand the occurrence of the confine¬ 
ment/deconfinement phase transition at finite temper¬ 
ature. Our result also confirms the well-known fact 
that the confinement/deconfinement phase transition sig¬ 
naled by the Polyakov loop average is associated to 
the center symmetry restoration/breaking and the in¬ 
finity/finiteness of the free energy for a single quark. 
Our derivation of the transition gives also the micro¬ 
scopic mechanism for confinement/deconfinement. In 


this sense, the gluonic mass at finite temperature is more 
important and powerful than the zero temperature, since 
it directly yields quark confinement. 

An important point missing in our analytical study is 
the lack of the analytical derivation of the gluonic mass 
in the same framework. At zero temperature, such a cal¬ 
culation has been given in the previous work [l3l| . More 
detailed investigations on the gluonic mass at finite tem¬ 
perature will be given in a subsequent paper [^ . 
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Appendix A: Integration over the fields 


The total Lagrangian except for the FP ghost term is 


-SfyM + -S^red + -2m — ~ 






(A.l) 


where we have defined 


+ (A.2) 


Here we have added the term to see the dif¬ 

ference between the massless and massive cases. In what 
follows, we consider the case of [P] =0, which means 
that is diagonal in gt, v. Then the total Lagrangian 
except for the ghost term and the higher-order terms in 
3T is cast into the quadratic form: 


- a:''^ 


1 

2 


c^A 









(A.3) 
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The integration over the field and jV is performed by 
the Gaussian integration and leads to the determinant: 

={det[{^p[r]9p[r])^^ + 

X det[(^'’['r]^p[y])^^ + (A.4) 

where the determinant over the Lorentz indices /i, is 
calculated. In the absence of mass M = 0, the deter¬ 
minant is independent of a. This is not the case in the 
presence of mass M ^ 0. The reduction condition is 
achieved at a = 0. Hence, we obtain the {D—1) massive 
vector modes and one massless scalar mode. 

The new variables have independent degrees of free¬ 
dom which are less than those expected from their ap¬ 
pearance. For instance, the field {A = 1, ...,dimG = 
— 1) has the independent components (a = 
1,..., dim(G/H)). In oder to integrate out the field vari¬ 
ables according to the path-integral formulation, we must 
correctly specify independent degrees of freedom for the 
respective field variable. This was achieved by adopting 
an appropriate reference of frame for the target space of 
the field [l^ for SU{2), and section 7.3 and Appendix 
J in for SU{3). Consequently, we can write the La- 
grangian in terms of the independent degrees of freedom 
as follows. 

Thus the SU{3) Yang-Mills Lagrangian density reads 

^YM = - 

+ 0(a3), (A.5) 

where we have defined 

Qt ■■= + 2ig^^,[r]{Tcr\ (A.6) 

with 

K<^b ^ _ gG^pgGP'^ 

+ [2gG^pdP + dPigG^p)]f‘^^. (A.7) 

For G = SU{3) in the m 2 Lximal option, is given in 
the usual Gell-Mann basis by 

j, fee {3,8}, a, &, c,de {1,2,4,5,6, 7}. (A.8) 

For G = SU{3) in the minimal option, is given in 
the usual Gell-Mann basis by 



j,fce {1,2,3,81, a,6,c,de {4,5,6,7}. (A.9) 

For G = SU{2), in particular, we have 

ia,b,c,d€ {1,2}), 
(A.IO) 


with 


_ ^_QPg^ ^ gGPgGpjS’^” + [2gGpdP + 9'’(gGp)]e“^ 

(A.ll) 


where we have used g“'=3e'=b3 _ ^ac^cb _ 

We introduce the complex-valued field $=*= defined by 


:= ± G C, 

y/2 


(A.12) 


and the inverse relation is given by 
4)^ = -I-$~), -I- i4)~). (A.13) 

V 2 V 2 


The term quadratic in 4>“ can be diagonalized by the 
complex variable: 


^^[-{S!p[y]9p[y])^^ + 

=^+[-{dp-igGpf+M^]^- 
+ ^-[-[dp + igGpf+ M^]^+. (A.14) 

Here we have used 

paqa ^ pagabQb ^ p+Q- p-Q+^ 

pa^abgb^ ■(p+Q-_p-Q+)^ (A.15) 

Thus we obtain 

det[-(^'’[r]^p[r])^® + 

= [-Dl[G]+M^][-Dl[G]+M% (A.16) 

where 

Dp[G] :=d- igGp, Dp[G] := 5 + igGp. (A.17) 

This yields 

Trln[-(^'’[r]^p[r])^^ + 

= lndet[-(^^[r]^p[r])^® + 

=Trln[-£>2[G] -f M^] + Trln[-£»2[G] + m"^]. (A.18) 
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Appendix B: Power series expansion and the coefficients 

The effective potential f4ff(v33, at i^s = 0 is given by 

Peff(¥>3,0)/r^=p-l) [f^(v33) + 2F^(^)] - [Fo(¥. 3) + 2Fo (^) 


(B.l) 


where 


d^-^p 


^m(¥’ 3) + 2-Fyf ((/33/2) — y ^27 ^)-D-i i¥’3)+2/^(p , p ‘ i / 2 )\ 


poo 

=Cd / dp p^-‘^[fj^{p^,(p3) + 2fj^{p‘^,P3/2)] 

Jo 

=Cd r dp p^-2{ln[l + e-2Vp=+^' - 2e-V^^^ cos(¥>3)] + 21n[l + _ 2e-v'p"+^" cos(¥>3/2)]}. 

Jo 


(B.2) 


The cosine is expanded into the power series in p 3 about ps = 47r/3 at which L = 0: defining a := ps — 47r/3 


cos ^ = - t - iCSiT + -tsrC S-“7=0^ - + O (cr^) . 

2 4 16 32^3 768 ^ ' 


2 


Then fj^{p^, P 3 ) is expanded into the power series in cr: 


ImIp , Ps) =cptip) + Cj^{p)a + Cj^{p)a + c%(p)a^ + Cj^{p)a + 0{a^), 


with the coefficients: 


c^(:p) =ln(l + e-^>’+e-2^>’) > 0, 


4(p) ■= - 


•\/3i 


e r 


1 + e + e 


< 0 , 


2 /ox e ^*>(1 + 4e + e 
Cm(p) = - < 0. 


4(p) = 


4(p) = 


2(1 + e~^p + e“2ep)2 
e-ep(l _ g-^P _ - e~^^p + e~^^p) 

2\/3(l + e~^p + e~^^p)^ ’ 

e-^p (1 + 12e-"*> - 12e-2'^p - 56e-^^p - 12e-^^p + 12e-^^p + e"®"'’) 

24 (1 + e“*^p + e~^'^p)^ 


where 


ip '■= yp^ + M2. 


The expansion of the integrand fj^{p'^, ps) + 2fj^{p‘^, p 3 /2) into the power series of a is given by 


(B.3) 


(B.4) 


(B.5) 


(B.6) 


fiyif, Pi) + V 33 / 2 ) (p) + h'^^J{p)a + h'-t’ {p)a^ + {p)a^ + h''y{p)a* + 0{a^), (B.7) 


,(i) 


(2)/ox 2 I I 1 ,( 4 ) / o \....4 


with the coefficients: 

x(0)/^3 -qJO) 


=^4 (p) = 31n(l + e + e > 0, 

h^V(p) =-c^^Mp) = _3 e~^"(l + 4e-^"+e~"^") 
V'T'l 4 (1+ e-^"p+e-2«p)2 


< 0 




73 e-"'7l - e-"" - 6e-2^'p - + e"'*"'’) 


(1 + e ’^p + e 2ep)3 


( 4 )sox _9 ( 4 )sox, _ 3 e-"r (1 + 12e-"r - 12e-^^p - 56e-^^p - 12e-^^p + 12e-5^'*> + e-6"»>) 




(1 + e + e ^^p) 


(B.8) 
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